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Abstract: In previous work we have shown that the equivariant index of multi-centered 
M = 2 black holes localizes on coUinear configurations along a fixed axis. Here we provide a 
general algorithm for enumerating such collinear configurations and computing their contri- 
bution to the index. We apply this machinery to the case of black holes described by quiver 
quantum mechanics, and give a systematic prescription - the Coulomb branch formula - 
for computing the cohomology of the moduli space of quiver representations. For quivers 
without oriented loops, the Coulomb branch formula is shown to agree with the Higgs branch 
formula based on Reineke's result for stack invariants, even when the dimension vector is not 
primitive. For quivers with oriented loops, the Coulomb branch formula parametrizes the 
Poincare polynomial of the quiver moduli space in terms of single-centered (or pure-Higgs) 
BPS invariants, which are conjecturally independent of the stability condition (i.e. the choice 
of Fayet-Iliopoulos parameters) and angular-momentum free. To facilitate further investiga- 
tion we provide a Mathematica package "CoulombHiggs .m" implementing the Coulomb 
and Higgs branch formulae. 
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1. Introduction 

BPS states in A/" = 2 supersymmetric string vacua offer a rich playground for exploring the 
microscopic properties of black holes in quantum gravity. Indeed, such states have dual de- 
scriptions as black hole solutions in supergravity at strong coupling, and brane configurations 
at weak coupling. The Witten index (more precisely, the second helicity supertrace[l, 2]) 
remains unchanged upon varying the string coupling. Computing the index ^2(7) in both 



- 1 - 



regimes as function of the electromagnetic charges 7 can provide non-trivial tests of the 
equivalence between the microscopic and gravitational descriptions expected in any theory 
of quantum gravity [3]. 

This comparison however is complicated by the fact that on the macroscopic side, con- 
tributions to the index ^(7) originate not only from single centered black holes with charge 
7, but also from multi-centered black holes with constituents carrying charges such that 
7 = ^ttj [4, 5, 6, 7]. The total index carried by single and multi-centered black holes is the 
quantity that should be compared with the index on the microscopic side. In contrast, the 
gravitational path integral in the near horizon geometry of each black hole gives information 
only about the index fis(a) associated with single centered black holes[8, 9]. Thus it is neces- 
sary to express the contribution from multi-centered black holes in terms of the index Qs{<^i) 
associated to each center. For this it suffices to compute the index^ fi'couiomb {q}) of 
the supersymmetric quantum mechanics of n centers carrying charges {ai} interacting by 
Coulomb and Lorentz forces (and other forces related by supersymmetry) . This 'Coulomb 
index' was computed in [10, 11], leading to a general prescription - the Coulomb branch 
formula - for expressing the total index in terms of the single-centered BPS invariants 
flsic^i)- Unlike these single-centered BPS invariants, the Coulomb index (7couiomb({ai}, {q}) 
depends on the moduli at infinity - indeed, this dependence is responsible for jumps of 
the total index ^2(7) across walls of marginal stability, providing a physically transparent 
derivation of the wall-crossing formulae known in the mathematical literature [12, 10, 13]. 

An important technical tool in [10, 11] was to consider a 'refined' version of these indices, 
which computes the trace Tr '(— l)^?/^"^-', where J3 is the angular momentum generator with 
respect to a fixed z-axis, and the prime denotes the removal of fermion zero mode contribu- 
tions before taking the trace. Unlike the helicity supertrace Tr'(— 1)^, the refined index is 
not a protected quantity in fuU-fiedged string theory and cannot be directly compared with 
the microscopic results. This generalization was nonetheless necessary, as one could use lo- 
calization with respect to rotations along the z axis to compute the refined Coulomb index 
fl'Couiomb({ai}, {cj}; y) for any number of centers n, and then take the limit y — )• 1 at the end 
to recover the result for Tr '(— 1)^. Although the description of the fixed points of this action 
is straightforward, their enumeration requires solving a set of n — 1 algebraic equations in 
n — 1 real variables, which quickly becomes unpractical as the number of centers increases. 
The first goal of this paper is to remove this bottleneck and give a completely algorithmic 
way of computing the Coulomb index (^couiomb {ci}',y), and hence the total refined in- 
dex n(7;y).^ Using the connection between multi-centered black hole quantum mechanics 
and quiver quantum mechanics described in [14] this leads to an explicit expression for the 
Poincare polynomial of quiver moduli spaces. The second goal is to establish the equivalence 

^The Fayet-Iliopoulos (FI) parameters {c^} depend on and on the values of the moduh at infinity, and 
will be treated as independent real parameters subject to the constraint = 0. 

^In [10, 11, 14] we used a subscript lof to denote a refined index. In this paper we shall drop this 
subscript to avoid cluttering up the formula:, but it should be understood that whenever the index carries 
the argument y, it corresponds to a refined index. 
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of this Coulomb branch formula with Reineke's formula[15] and its generalizations [16, 12] for 
the Poincare polynomial of quiver moduli spaces for general quivers with no oriented loop. 



Coulomb index from localization 

Before explaining our new algorithm, let us briefly review the prescription of [11] for com- 
puting the refined Coulomb index gcouiomhiic^i}, {ci}',y)- Let T be the charge lattice of 
electromagnetic charges, equipped with the Dirac-Schwinger-Zwanziger (DSZ) symplectic 
product (■,■) G Z. Consider a multi-centered black hole configuration where each center 
carries charge ai {i = 1. . .n), with DSZ products atj = {ai,aj) between the charges. For 
fixed values of the moduli at infinity, encoded in the FI parameters q, n-centered configu- 
rations are parametrized by a 2r;, — 2-dimensional phase space Ain{{c(ij}, {ci}) [4, 5]. The 
coordinates of A^„ are the locations of n centers rj up to overall translations, subject to 
n — 1 constraints^ 

Vz = l...n, V^L=c.. (1.1) 

The space Ain admits a symplectic form [5, 17, 18], such that the action of S0{3) rotations 
on is generated by the moment map 

equal to the angular momentum carried by the configuration. We denote by 

5'Couiomb({ai,---an};{ci,---c„};y) = Tr'(-?/)^'^3 (1.3) 

the refined index associated to this multi-centered configuration assuming that all the centers 
are distinguishihle from each other, and that each center carries no intrinsic degeneracy. 
Mathematically, this is the equivariant Dirac index of the symplectic space A^„({ajj}, {q}) 
[11, 19]. We refer to (1.3) as the Coulomb index of multi-centered black holes with charges 
{ai}. This in turn can be used to compute the index associated with a general multi-centered 
black hole configurations in terms of indices VL^{ai) carried by individual centers following 
the procedure described in [11, 14]. 

When there exists an ordering of the charges {aj} such that i < j ii and only if aij > 
0, and for generic values of the parameters q away from walls of marginal stability, the 
symplectic space Ain is compact, and the Coulomb index can be computed by localization 
with respect to rotations around a fixed axis, using the Atiyah-Bott Lefschetz fixed point 



■^The constraint (1.1) for i ~ n follows from the sum of others using the fact that X^ILi '^i ^ ^"-"^ 
a general multi-centered black hole system there are additional constraints besides (1.1) coming from the 
requirement of the regularity of the metric and other fields. However when the central charges of the 
constituents nearly align, which is the the limit in which the quiver quantum mechanics becomes a good 
description, these additional constraints arc expected to be satisfied automatically [14]. Throughout this 
paper we shall be working in this limit. 
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formula [10, 11]. The configurations which stay invariant under such rotations are colhnear 
configurations, corresponding to solutions of (1.1) lying along the z axis. The result is 
expressed as 

^?Coulomb({«l,---«n};{Cl,---C„};y) = (-1)"-1+E,<,-. (^_^-l)-+l ^ ± ^E.<,-. sign(.,-.,) ^ 

extrema 

(1.4) 

where the sum runs over solutions to the equations 

n 

^ =Ci, forl<z<n-l, zi = 0. (1.5) 

The = condition fixes the translational zero-mode. The sign ± in (1.4) is given by the 
sign of the Hessian ^QiiS^W j dzidz^ of the superpotential 

n 

= - ^^Sn(2:j - Zi) log \zi - Zj\ - ^^CiZi, (1.6) 

i<j 1=1 

whose critical points reproduce the conditions (1.5). 

For charges {tti} such that no such ordering exists, the space Ain may be non-compact, 
due to the possibility of 'scaling solutions', i.e. a subset of the centers approaching each other 
at arbitrary short distances [20, 7]. In that case, we continue to define the Coulomb index 
fi'Couiomb({tti5 ■ ■ ■ ttn}; {^1, ■ ■ ■ c„}; y) hj the localization formula (1.4), although the result can 
no longer be interpreted directly as the refined index associated to the multi-centered black 
hole configuration (in particular, it may not be a symmetric Laurent polynomial). Never- 
theless it can be used to construct such a refined index following the procedure described 
in [11, 14] and reviewed in §3.1. When the FI parameters Cj sit on a wall of marginal sta- 
bility, the space is also non-compact due to the possibility of separating the centers 
into two or more clusters at arbitrarily large distances, and we do not assign a value to 
fi'Couiomb({ai, ■ ■ ■ an}; {ci, ■ ■ ■ c„}; y) in such cases. 

A new algorithm for computing the Coulomb index 

Except in very special cases, it is usually impossible to find all solutions of (1.5) explicitly. 
This is also unnecessary, since the contribution of a given solution of (1.5) to the total 
index (1.4) depends only on the ordering of the centers, via the angular momentum J3 = 
i Si<j '^ij sig^(^j ~ ^i) along the 2;-axis and the sign of the Hessian W". For a small number 
of centers, it is possible to find approximate solutions numerically, and determine both J3 
and the sign of W", however this becomes quickly unpractical as the number of centers grows. 
Moreover, the brute force enumeration of solutions of (1.5) does not take into account the 
fact that there can be cancellations between different solutions with the same ordering. To 
exploit this fact, it is useful to associate a permutation a to each solution to (1.5), such 
that z < j iff Zo-(j) < z„(^jy Defining cti = a„(^i), Xi = Za-[i) and q = Co-(i), solutions of (1.5) 
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correspond to critical points of 

n n 

W{{xi}) = - ^aij\og\xi-Xj\-^CiXi, ^Q = (1.7) 

i<j i=l i=l 

in the physical region 

= Xi < X2 < X3 ■ ■ ■ < Xn ■ (1.8) 

Reorganizing (1.4) as a sum over all permutations a of 1, 2, ... n, we obtain 
^?Couiomb({«i, ■ ■ ■ an}; {ci, " " " Cn}; v) = (-l)"-i+S»<. (y - y-')-''+' s{a) yS«<.-'^w^U) , 

(1.9) 

where s{a) is the sum of the sign of the Hessian of (1.7) over each critical point in the 
physical region. In particular, s{a) is insensitive to pairs of solutions of (1.5) with the same 
ordering, which may appear under small perturbations of the parameters aij and q, as long 
as we stay away from walls of marginal stability in the space of FI parameters {q} and 
away from certain 'scaling walls' in the space of DSZ products {aij} described in more detail 
below. 

The first aim of this paper is to develop an explicit algorithm for computing s{a) and 
hence the Coulomb index (1.9) for generic DSZ products {aij} and FI parameters {cj}. 
This is achieved in §2, where we prove an inductive formula (2.20) for the indexed number 
s(cr) = F{{ai, ■ ■ ■ an}, {ci, ■ ■ ■ c^}) of critical points of the superpotential (1.7), by exploiting 
its robustness under changes of the DSZ products aij. For quivers without oriented loops 
and generic products, we arrive at the following completely explicit result: 



^Couiomb({ai,---a„};{ci,---c„};2/) =(-1)" i+J:»<. (y - y 1) 

n-1 / k \ 

5^n® h'^w.-e^^+i) 5^c.(i) (-i)^^-^'®(-"''«-"'+^')2/^ 



;i.io) 



fc=i 



where 6 is the Heaviside function and the sum runs over all permutations cr of 1, 2, ■ ■ ■ n. If 
some of the aj/s vanish then we need to deform them away from zero such that the deformed 
quiver continues to satisfy the no loop condition, compute the result using (1.10) and then 
take the limit where the deformations are taken back to zero. 

As described in (2.20), in the presence of loops the expression for s{a) picks up additional 
contributions AFa given in (2.32), which depend on the index F with fewer centers and 
another auxiliary quantity G{ai, ■ ■ ■ an)- The latter counts the (indexed) number^ of coUinear 
scaling solutions, i.e. solutions of (1.5) with Cj = which may arise when the total angular 
momentum ^ Yli<j ^ij vanishes. We find that this 'scaling index' can itself be computed 
inductively using (2.38). These formulae hold when the DSZ products otij are generic, but 
we show that even when this is not the case, all the relevant physical quantities can be 
computed in terms of limits of these formulae. 

''Due to the scaling symmetry, each solution arises as a one-parameter family, the number of which is 
counted by G. 
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Quiver quantum mechanics and pure-Higgs invariants 



While the refined index is not protected in full-fiedged string theory, it is protected in the 
context of A/" = 4 supersymmetric quiver quantum mechanics, which describes the dynamics 
of certain multi-black hole bound states around special loci in moduli space where the central 
charges of the constituents become nearly aligned [5, 7]. Thus, by considering a black hole 
whose dynamics in some region of the moduli space is described by a specific quiver quantum 
mechanics, we can use our general Coulomb branch formula for multi-black hole bound states 
to parametrize the refined index in the corresponding quiver quantum mechanics in terms of 
single-centered BPS invariants [14]. For quivers without oriented loops, the single-centered 
BPS invariants are trivial, and the Coulomb branch formula is completely explicit. We can 
then use this to establish the equivalence of the Coulomb branch formula for such quivers 
with explicit formulae for the cohomology of the moduli space of stable quiver representations 
known in the mathematical literature [15, 16, 21, 22]. 

Before explaining our results, let us briefiy review the relation between quiver quantum 
mechanics and multi-centered black holes [14]. A/" = 4 supersymmetric quiver quantum 
mechanics can be obtained by dimensionally reducing an A/" = 1 supersymmetric gauge 
theory in 3-1-1 dimensions^ - containing vector multiplets in the adjoint representation of the 
gauge group G = Y[i=i K^i^i) ^^"^ chiral multiplets in bi-fundamental representations of 
U{Ni) X U{Nj) - down to 0-1-1 dimensions. The scalars coming from the chiral multiplets 
are called the Higgs variables and those from the vector multiplets are called the Coulomb 
variables. The vacuum moduli space of the Higgs variables at zero values of the Coulomb 
variables is equivalent to the moduli space Ai of stable quiver representations (in short quiver 
moduli space), where the stability condition is determined by the FI parameters (i, . . .(k 
for each U{Ni) factor. BPS states are in one-to-one correspondence with cohomology classes 
in H* {A4 ; Z) , and the angular momentum is identified with the component J3 = {p — d)/2 
of the Lefschetz SU{2) action on the total cohomology H*{Jli; Z). Thus the refined index is 
given by the 'Poincare-Laurent polynomial'^ 



where bp{Ai) are the Betti numbers and d is the complex dimension of Ai. But the same 
spectrum may also be calculated by first integrating out all the Higgs variables and consid- 
ering the effective theory for the Coulomb variables. The latter turns out to be given by 

^As an aside it should be noted that A/" = 4 supersymmetric quiver quantum mechanics is also useful in 
computing the spectrum of BPS states in Af = 2 supersymmetric gauge theories in 3-1-1 dimensions [23, 19]. 
Indeed, the BPS spectrum of many gauge theories can be understood in the language of quiver representations 
[24]. For example in the context of A/" = 2 supersymmetric pure SU(2) gauge theory, the role of single centered 
'black holes' is played by the monopole and the dyon in Seiberg-Witten theory, which are stable for all values 
of the moduli and whose bound states generate the complete BPS spectrum. 

^Wc use this terminology since Q{A4; y) differs from the usual Poincarc polynomial bp-y^ by a y — > — y 
transformation and a multiplicative factor of (— y)^^. 
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the same quantum mechanical system as that of multi-centered black holes in A/" = 2 su- 
persymmetric string theory, upon identifying the charge vector 7 with the dimension vector 
{Ni, . . . Nk) [5], and the DSZ product 7,^ between the basis vectors 7^ = (0, . . . , 1, . . . , 0), 
where the only non-vanishing entry occurs in position z, with the number of arrows from 
the 2-th to the j-th node of the quiver. The Coulomb branch formula of [11] can thus be 
used to express the Poincare-Laurent polynomial (1.11) in terms of certain 'single centered 
BPS invariants' r2s(a), which are (conjecturally) independent of the FI parameters and of 
the fugacity parameter y [14]. In simplest cases, ^^3(7) enters just as an additive constant 
in Q{Ai;y) = ^(7), and can be identified as the Lefschetz singlet contribution to the total 
cohomology iJ*(A^;Z) [25, 26, 14, 27]. In general however, the single centered BPS invari- 
ants fis(a) enter in Q{A4.\y) in a more complicated fashion. It is a very interesting open 
problem to identify the corresponding 'absolutely stable' classes in H*{Ai; Z). 

For quivers without oriented loops, the only non- vanishing single centered BPS invariants 
are those associated to the basis vectors 7^, and for such vectors f2s(7i) takes the value 1. 
In that case, the Coulomb branch formula gives a completely explicit result for Q{Ai;y), 
which can be compared to known results in the mathematical literature. In particular, for 
primitive dimension vector 7 (i.e. such that all Ni are coprime), Reineke's formula [15] gives 
another completely explicit result for Q{Ai;y). For Abelian quivers (i.e such that all A', 
are one), the Coulomb branch formula equates the Poincare-Laurent polynomial Q{A4.]y) 
with the Coulomb index 5'couiomb({7i, ■ ■ ■ In}', {Ci; ' ' ' Cn}; y) given in (1.10). In §4.2 we show 
the equivalence of the Reineke's formula and Coulomb branch results for Abelian quivers, 
generalizing previous arguments given in the context of wall-crossing [10, 28]. For non- 
Abelian quivers with primitive dimension vector, the equivalence of the Reineke's formula 
and Coulomb branch results can be reduced to the Abelian case, by using the Abelianization 
property satisfied by Reineke's formula [28, 21] (see (4.18) below). 

For non- Abelian quivers with non-primitive vector, Reineke's formula no longer computes 
the Poincare-Laurent polynomial Q{Ai;y) of the quiver moduli space Ai (which is singular 
due to marginal bound states), but rather the 'stack invariant' Gniggs- It is conjectured 
in the mathematical literature that a bone-fide Poincare-Laurent polynomial Q{Ai; y) may 
be reconstructed from the stack invariants Gniggs [16, 12] (see Eq. (4.1) below), but it is 
unclear in general how to construct a smooth moduli space Ai whose cohomology would 
agree with Gniggs- At any rate, using a further property of Abelian stack invariants (^Higgs 
established in [22] (see (4.22) below), we prove that the Coulomb branch formula agrees with 
the Poincare-Laurent polynomial Q{A4;y) computed from Reineke's formula for the stack 
invariants. The result can be written as 




/i(m) y-y 



1 




(1.12) 



-7- 



where /i(m) is the Mobius function, {(^7j) ^ } denotes that we have kj nodes each carrying 
charge £7-,- for £ > 1 and I < j < n, and it is understood that in computing (^couiomb whenever 
some aij vanishes we need to deform it away from that value to produce an Abehan quiver 
without loop, and then use (1.10) for computing it. 

For brevity we shall henceforth refer to Reineke's formula and its generalizations as the 
Higgs branch formulae, since they compute directly the cohomology of the Higgs branch 
moduli space. For quivers with loops, the Higgs branch formula will refer to the result of 
any computation of the cohomology of the Higgs branch moduli space, although no general 
formula is available in such cases. 

A Mathematica package for quiver invariants 

It should be clear from the above that the Coulomb branch formula produces a parametriza- 
tion of the Poincare-Laurent polynomial of any quiver in terms of single-centered BPS invari- 
ants in a completely combinatoric way. However, even for moderately complicated quivers the 
necessary computations quickly become tiresome, and are best implemented on a computer. 
We have implemented the Coulomb branch formula as well as Reineke's formula for stack 
invariants (and many related other routines) in a mathematica package "CoulombHiggs" 
available from arXiv and described in Appendix A, which we hope will facilitate studies of 
single-centered BPS invariants. This package has been successfully tested on the examples 
investigated in [14] and many more. 
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2. A formula for the Coulomb index of multi-centered black holes 

In this section, we establish a recursive algorithm for computing the Coulomb index (?couiomb({«i, ■ ■ ■ 

{ci, ■ ■ ■ c„,}; y) for general charge configurations a,, away from the walls of marginal sta- 
bility in the space of FI parameters Cj. We start in §2.1 with charge configurations for which 
there exists a possible ordering of the aj's such that 

aij > for i<j, (2.1) 

and with all DSZ products a^j non-zero. It follows from the discussion in §1 that this 
corresponds to an Abelian quiver without any oriented loop. For such systems we obtain the 
simple formula (2.9) for the indexed number s{a) that enters (1.9). In §2.2 we show that for 
the purpose of computing the total Coulomb index, the same result can be used even when 
some of the a^j's vanish. In §2.3, we turn to general multi-centered black hole configurations 
for which the charges do not admit an ordering satisfying (2.1), and establish an inductive 
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formulae for computing s{a). This is given by F{{ai, ■ ■ ■ an}, {ci, ■ ■ ■ c„}) in (2.20) with 
tti = ao-(j). During this analysis we also derive a similar formula for the coefficient s(cr) for 
coUinear scaling solutions[20, 7] for which all the FI parameters vanish and the a^'s satisfy 
S«<j ~ 0. The corresponding inductive formula for s(cr), called G(ao-(i), ■ ■ ■ aain)), 

is given in (2.38). Combining these results and summing over all permutations yields a 
general algorithm for 5'couiomb({«i, ■ ■ ■ «n}; {ci, ■ ■ ■ c„}; y) with non vanishing DSZ products, 
as summarized in (2.33). 

2.1 Abelian quivers with no oriented loops and all aij ^ 

We start by considering the case of charge configurations which admit an ordering such that 
aij > for i < j. A special case arises if all ai are positive linear combinations of two charge 
vectors 71,72 with (71,72) > 0, as is the case for multi-centered configurations relevant for 
wall-crossing [10]. 
Defining 

n i 

Hi = Xi+i -Xi, di = - ^ Ci = ^Ci, for 1 < z < n - 1 , (2.2) 

j=i+l j=l 

we can express (1.7) as 

i<i \k=i / i=l 



dW ^ 



This gives 

^ = -5^ 5^ a,,-^ + 4 forl<A;<n-l. (2.4) 

Each yi takes value from to 00, but by an appropriate coordinate transformation Ui = f{yi) 
we can bring the range to < Uj < 1 for each i. We can regard the space spanned by the 
Uj's a unit box. Our goal is to examine the condition under which W has an extremum 
with respect to the ?/j's in the interior of this box, i.e. there is a solution to the equation 
dW/dyk = for every k. 

To this aim, let us now consider the following deformation of the cij/s: 

ai,i+i ai,i+i Vz, aij Xoij for |z - j| > 2 , < A < 1 . (2.5) 

In the limit A — 0, only cij j+i's remain non-zero and (2.4) takes the simple form 

^— = + 4. (2.6) 

oyk yk 

Thus the set of equations dW/ dy^ = has a solution in the range < j/^ < 00 if and only if 

sign(afc,fc+i) = sign(4) for 1 < < n - 1 . (2.7) 



Furthermore the sign of the Hessian of W at this solution is easily determined to be 

n-l 

n sign(afc,fc+i) . (2.8) 

k=l 

These results can be summarized by saying that for quivers without loops the coefficient 
s{a) of y^i<j°'^3 associated with a given permutation {ai, ■ ■ ■ a„} is given, for A = 0, by 

n-l 

Fo{{au ■ ■ ■ ■ ■ ■ = n 4)(-l)^^-' e(-5,,,+o ^ (2.9) 

fc=i 

where G(x) is the Heaviside function. 

Now consider changing A from to 1. During this deformation new extrema can appear 
and disappear in pairs in the interior of the box, but since they are weighted by the sign of 
the Hessian, the result is unaffected. The other possibility is that new extrema can emerge 
from (or disappear into) the boundary of the box. This happens when a subset of the |/i's 
approach and/or infinity. Since dW/dyk approaches the constant dk a.syk ^ oo irrespective 
of the values of the other ?/j's, it is clear that away from the marginal stability walls dk = 0, 
none of the |/j's can approach infinity. Thus the only possible boundary component where 
extrema of W can appear or disappear is where a subset of the ?/j's vanish. 

Let us suppose that as A is increased from to 1, such a phenomenon takes place at 
some value A = Ac- If a new extremum appears at Ac then for A slightly above Ac, there 
will be an extremum of W where a subset of the yi's are small, corresponding to a subset 
of the centers being close to each other. If on the other hand an extremum disappears as A 
approaches Ac from below, then such a configuration exists for A slightly below Ac. These 
correspond to onset or disappearance of scaling solutions[20, 7], with A = Ac being the point 
at which the scaling solution becomes collinear. 

Now if the quiver corresponding to the original a^^s had no oriented loop then neither 
does the quiver associated with the deformed ajj's, since the signs of all the ajj's are preserved 
under the deformation (2.5). This implies that the deformed quiver cannot have a scaling 
configuration, and hence, as we deform A from to 1, no extremum of W can emerge from 
or disappear into the boundary. Thus (2.9) gives the correct contribution to s{a) even at 
A = 1. Using (1.9) we now get 

<7couiomb({«i, ■ ■ ■ «n}; {ci, ■ ■ ■ y) = (- l)"-i+S»<. {y - y-')--^' 

n-l / k \ 

5^ n ® J2 (-1)^^-1 ^E.<,a.W.{i) 

a k=l \ i=l / 

2.2 Abelian quivers with no oriented loops but some aij = 

We now turn to the case of quiver without oriented loops, but for which some of the aj/s 
vanish. For this we first deform all the vanishing a^j's to non-zero values in such a way that 
the deformed quiver does not have any oriented loop. To see that this is always possible, let 



(2.10) 
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us carry out the deformation one link at a time. We begin with the original quiver without 
oriented loop and make one of the vanishing ai/s non-zero. If this leads to a quiver with an 
oriented loop, then there exists some component C of the original quiver, which, together 
with the new link, gives rise to a quiver with oriented loop. If so let us flip the sign of aij of 
the deformed link. In this case C together with the added link no longer forms an oriented 
loop. Suppose there were another component C" of the original quiver, which combined with 
the new link would now form an oriented loop. Then C + C' would form an oriented loop in 
the original quiver, whih contradicts our assumption. Thus, by choosing the sign of a^j of 
the deformed link we can ensure that the new quiver also does not have any oriented loop. 
We can now repeat the argument and show that all the vanishing a^j's can be made non-zero 
and for appropriate choice of sign of the deformed aj/s the new quiver does not have any 
oriented loop. Thus we can compute its index by our earlier formula (2.10). 

We shall now argue that the index of the original quiver can be obtained by taking the 
limit of the index of the deformed quiver in which the deformation parameters go to zero. 
For this we shall work with the total index (1.9) rather than a given permutation. We shall 
use the original charges and use their locations Zi - related to the Xi by Xi = Zo-(i) - as the 
independent variables. The ^j's satisfy (1.5). 

Let us now consider the effect of taking apg to for some specific p, q. For any extremum 
of W at which the locations of Op and ag remain at finite separation, this limit has no drastic 
effect and the contribution to the index from such extrema at apg = is the same as what 
we get by taking the apg — )■ limit. Thus we only have to examine the fate of the critical 
points for which the locations of ap and ag approach each other in the apg — )■ limit, as 
generically such critical points will disappear for apg = 0. For such solutions we can replace 
Zg by Zp in (1.5) except in the apg/\zp — Zq\ terms, and express (1.5) as 



- E 



Oiik oipk + agk 



_^ ^k] ^k 



— Cfc =0 ioT 1 < k < n, k p, q 



(2.11^ 

E(-^ip ^qp |-, \ ^iq ^pq p. 

I i" ~ I i " ~ ' ~ / > I i" ~ I T ~ ^9 ~ U . 

\7i — 7 \ r — 7. ^ \7- — r. \7. — 7 \ 



i—i ^p\ \^q ^p\ i—i ^p\ \^q 

By adding and subtracting the last two equations we get 



E 



Olik Olpk + Olqk n t ^ ^ 1 ^ 1 u 

Cfe =U lOT 1 < k < n, k p, q 



\Zp Zk\ 

n 



y a,p + a,g _ Q ^ y a..p a,g ^ ^ apg _ q _ 

^ ^ y . 7 \ ^ 7 ■ 7 \ \ 7 7 \ 



i^i i=i \^<1 



(2.12) 



The first set of equations and the second equation together correspond to the equilibrium 
configuration of n — 1 charges in which the charges ap and ag have merged to form a charge 
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Op + aq and the corresponding FI parameters have been added. The last equation can be 
interpreted as an equation for Zq — Zp. The existence of a solution to this equation requires 



sign 



Cq 



^ I Zp I 



signttqp. (2.13) 



When this condition is satisfied then Zq — Zp is of order apq for small apq. On the other hand 
when apq = 0, the last equation in (2.12) generically has no solution since the left hand side 
of the equation becomes independent of Zq and all the Zj's for i ^ q are already fixed by the 
other equations. This shows that the critical points of W associated with solutions to (2.12) 
disappear at apq = 0. This could give rise to a discontinuity in the index at apq = 0. The 
important point to note however is that even if (2.13) is satisfied, the solutions to the last 
equation in (2.12) always occur in pairs, related by a reversal of the sign of Zp — Zq. The 
exponent of y in (1.4) remains unchanged under this exchange in the apq — )■ limit since 
this exchange only flips the sign of the coefficient of apq in the exponent. Finally it is easy 
to check that s{a) changes sign under this exchange. Thus the contribution from this pair 
of solutions cancel and we get a smooth apq — > limit. Repeating this analysis for the other 
apq^s we see that the index associated with the original quiver can be obtained as the limit 
of the index associated with the deformed quiver. 

Note that the above argument breaks down if the solution to the equations in the first 
two lines of (2.12) automatically satisfy 



Cq Cn 



y o^^P "^:' =o. (2.14) 



! = 1 



This happens for example when ap and aq are parallel so that aipja^q = \ap\/\aq\ = Cp/cq. In 
this case at apq = there is a solution to (1.5) at Zp = Zg, obtained by solving the equations 
in the first two lines of (2.12). Now consider the case when apq is deformed away from 0. In 
this case in order to look for a solution to the last equation in (2.12) where Zp and Zq are 
close to each other we can no longer set Zp = Zq in the regular terms from the beginning, 
but must keep terms of order {zp — Zq) in the last equation. If we call this term A{zp — Zq) 
for some constant A then we can express this equation as 

A{zp-Zq) + 2 ^ =0. (2.15) 

This equation is no longer invariant under a change of sign of Zp — Zq, and in fact has a 
solution only for one sign of Zp — Zq irrespective of the sign of apq. In the apq — )■ limit this 
solution smoothly continues to the solution with Zp = Zq at apq = 0. Thus we again see that 
the apq limit is smooth, and agrees with the result for apq = 0. 

2.3 Generic Abelian quivers with all aij non-zero 

We shall now consider a generic multi-centered black hole configuration with all aij non-zero, 
but whose associated quiver may possess oriented loops. In this case we need to take into 
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account possible contributions from scaling solutions. Our goal in this section will be to 
compute (7couiomb for such configurations. 



2.3.1 An inductive formula for the index of collinear solutions 

We proceed as in §2.1 and consider the deformation (2.5). It is clear that at A = the 
contribution from a given permutation a will be given by (2.9). Thus we need to investigate 
the total change in this contribution as A changes from to 1. As discussed in §2.1 these 
changes could come from values of A at which a set A of neighbouring centers come close to 
each other. This can happen if the total angular momentum carried by this set of centers 
vanish, 

E "^.=0' (2.16) 

where a'^j denotes the deformed aij. To see this, note that in the limit where all Ui for 
i,i + 1 G A approach zero, the part of the superpotential (2.3) involving the ?/j's becomes 
quasi-homogeneous , 

W{{\y,})^W{{y,})- J2 %log^- (2-17) 

i,j&A;i<j 

Differentiating with respect to A and using dW/dyi = implies (2.16). Since the set A must 
contain at least three elements there are at most [n — 2){n — l)/2 possible sets A, given 
by the {n — 2){n — l)/2 possible ways of choosing the beginning and the end of the set. 
Correspondingly there are at most (n — 2)(n — l)/2 possible values A^ of the deformation 
parameter A where such collinear scaling configurations can arise. Using (2.5), the condition 
(2.16) becomes a linear equation in A^i, 

Xa Y1 Y1 "M+i=0. (2.18) 

«jGA;i<j-2 i£A,i+leA 

The index can jump across A = A^ if A^ hes between and 1. This is so if and only if the 
left hand side of (2.18) has opposite signs at A^i = and 1, i.e. 



A,i+l£A / \i,j£A,i<j / 



, , ^ , (2.19) 

\i&A,i+l<£A 

If F({ai, ■ ■ ■ {ci, ■ ■ ■ Cn}) denotes the coefficient s(cr) of y^^<o associated with a given 
permutation, then we have 

F({q!i, ■ ■ ■ din}, {Ci, ■ ■ ■ Cn}) = i^o({«l, " " " "n}, {Cl, " " " C„}) + ^ AF4 , (2.20) 

A 

where Fq is given by (2.9), and AFa is the jump across the critical point A^. Our goal will 
be to compute the expression for AFa- 

Let us suppose that set A consists of the integers k, k + 1, ■ ■ ■ i. We shall examine the 
configuration close to the critical point by taking 

E ^rs=' (2-21) 



r.s 
k<r<s<i 
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for some small number e. We now define Zg via 

Xs = Xk + y Zg foTk<s<i, Zk = , ze = 1 , (2.22) 

and use xi, - ■ ■ Xk, xg+i, ■ ■ ■ x„, y, Zk+i, ■ ■ ■ -z^-i as independent variables. Then the relevant 
equations are given by extremizing 

k-i e 



W = - ^ a'ij\n{xj - Xi) -'^^a'ijn{xk + yzs- Xi) 

^,3 i=l s = k 

n e 

- ^ ^ a,i \n{xi -Xk- yzs) - ^ a,^ ln(i/(2r - Zs)) 

- ^ CiXi- \^^Cs\xk-y CsZs. (2.23) 



i<k or i>l + l 



^s=k / s=k+l 



We shall be examining an extremum of W for which y is small, of order e. In this case the 
extrema of W with respect to xi, - ■ ■ Xk, xe+i, ■ ■ ■ x„ can be obtained by extremizing 



fe-i / e \ 

Wi = - a[j \n{xj - Xi) - Yi^aiA ln( 

',i i=l \s=k / 



Xk Xi) 



i<Cj;i,j<k or >^+l 



- 5Z I 5Z "'^M - Xk) - ^ CiXi - j ^ c J Xfc. (2.24) 

i=l+l \s=k J i \s=k / 

i<k or i>e+l 

The existence of an extremum of Wi is equivalent to the existence of a coUinear configura- 
tion with n — £ + k centers with charges a[, ■ ■ ■ a'k^i, Yll=k ^s' ' ' ' FI parame- 
ters Ci, ■ ■ ■ Cfc_i, Yll=k Q+i, ■■■£„, situated at Xi, ■ ■ ■ Xk-i,Xk, x^+i, ■ ■ ■ Xn- Extremization of 
(2.23) with respect to the parameters Zs for k + 1 < s < i—1 can be obtained by extremizing 

W2 = - Yl <H^r-Zs). (2.25) 



s,r 
k<s<r<e 



The existence of an extremum of W2 is equivalent to the existence of a coUinear scaling 
configuration of i — k + 1 centers with charges a'f^, - ■ ■ a'^ and zero FI parameters, with the 
locations of the centers being at Zk = 0, Zk+i, ■ ■ ■ , ze-i and zi = 1. Finally y can be obtained 
by extremizing^ 



fc-l £ n i 



.^"'""^^ i=l s=k+l ^ * i=l+ls=k+l * s=k+l 

k<s<r<£ 



(2.26) 



''Since y is small, we have expanded the terms in W which arc non-singular in the y — > limit to first 
order in y. 
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Using (2.21) this gives 

fc-l I n £ e 

- E E a;.^-^ + E E - E = •) ■ (2-27) 

^ i=l s=k+l i=£+l s=k+l s=k+l 

For generic aij, a solution to this equation with positive y exists only for one choice of sign 
of e. We shall assume that we have taken the sign of e to be such that the solution exists. 
Let rjA denote a quantity which takes value 1 (—1) if the solution exists for A above (below) 
the critical value given in (2.18). 

Let us now compute the Hessian at this critical point. From (2.24)-(2.26) it is clear 
that the second derivative of W with respect to all the variables except y remain finite, and 
we have d'^W/dy'^ ~ Yl "•'^ ^srllP' ~ ^llP'- Since this is large for ?/ ~ e and all other 

fc<s<r<^ 

second derivatives of W remain finite, the full determinant will be given by the product 
of d'^W jdy^ and the determinant of the Hessian involving the rest of the variables. Fur- 
thermore we see from (2.23) that d'^W j dx^dzg goes to zero as ?/ — !■ 0. Thus the Hessian 
of W with respect to Xj's and Zg-s, factorizes into the product of the Hessian of W\ with 
respect to Xi, ■ • ■ x^, x^+i, ■ ■ ■ x„ and the Hessian of with respect to ■ ■ ■ Zi_x. In our 
notation the sign of the Hessian of W\ with respect to Xi, ■ ■ ■ x^, x^+i, ■ ■ ■ x„ is given by 
F({«'i, ■ ■ ■ Q!$,_i, a'fc H a^, 61^^^, ■ ■ ■ a^}, {ci, ■ ■ ■ Cfc_i, 4 H q, c^+i, ■ ■ ■ c„}). Let 

G(a'„ ■■■«;)= sign f det [d^A^'^S] (2-28) 

be the sign of the Hessian of W2 with respect to -Zfe+i, ■ ■ ■ zi_x when the corresponding scaling 
solution exists; otherwise we take GicJ^,, ■ ■ ■ a'^) = 0.^ Then we can write 

AFa = F{{a[, ■ ■ ■ a'f._^, a'^. -\ a^^^, ■ ■ ■ a^}, {ci, ■ ■ ■ Cfc_i, 4 H c^, 4+1, ■ ■ ■ c„}) 

i=k 

k<i<j<e 

where the last factor imposes the constraint (2.19). Now it follows from (2.5) that at A = 
\a + S\, 

e 

e= a'^^ = SX diij . (2.30) 

fe<r<s<f i<j-2 

Suppose the solution exists for 6X > 0. Then we have rjA = 1 and we see from (2.30) that 
sign(e)=sign(^\j=fe cijj). On the other hand if the solution exists for (5A < then we have 

i<j-2 

TjA = —1 and sign(e)=— sign(^ ij=fc aij). Thus in either case 

i<j-2 

I 

?7Asign(e) = sign(^ ^ oii^ . (2.31) 

i,j — k 
%<]~2 

*As usual if there are more than one solutions then we add their contributions. 
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Substituting this into (2.29) we get 
AFA = F{{a[,---ai_„ai + 



^+1, ■ ■ ■ {Ci, ■ ■ ■ Cfc_i, Cfc + ■ ■ ■ Q, ■ ■ ■ Cn}) 





( ' 


\ 




/ 




X (^(a'fc, ■■■«£) sign 


E 




e 














V 





E 

k<i<j<t 



\ 



J 



\ 



I 



(2.32) 

Note that a special case of (2.32) is /c = 1, £ = n in which case the F on the right hand 
side of this equation is -F(a'i + ■ ■ ■ a^; ci + ■ ■ ■ c^) = 1. If we can compute G{a.'f,, ■ ■ ■ «^), then 
we can use (2.20) and (2.32) to compute the function F recursively. Once we know how to 
compute F, the Coulomb index can be computed as 

5'Coulomb({ai, ■ ■ ■ «n}; {Cl, ■ ■ ■ Cn}; v) 



=a(n)}) y 



(2.33) 



where the sum runs over all permutations a. 

2.3.2 An inductive formula for the index of scaling collinear solutions 

We now turn to the computation of G{ai, ■ ■ ■ am), the indexed number of critical points of 
the superpotential 



W 



a 



0, 



(2.34) 



1 <z< j <m 



1 < z < J < Tn 



in the range Zj > Zi for j > i. This coincides with (2.25) under the identifications 
= {a'l^, ■ ■ ■ a'^} , and obvious redefinitions of Zi. The invariance of the super- 
potential (2.34) under both translation and rescaling of the Zi^s, must be 'gauge fixed' before 
counting critical points. These invariances were fixed by the conditions zi = 0, Zm = 1 in 
(2.22); however in order to compute G{ai, ■ ■ ■ am) inductively, it will be more convenient to 
choose a different gauge zi = 0, Zm-i = 1- 
Let us now consider the deformation 

m— 1 



Ciim &im fof i = 1, 2, ■ ■ ■ TTi — 1, a. 



m— 3,m— 1 



— )■ a 



m—3,m—l 



A)E 



cy.4. 



(2.35) 



SO that the deformed ai/s (called a'^j) continue to satisfy ^j<j a[j = 0. In the limit /i — )■ 0, we 
can treat the m-th center as a probe in the background of other centers. From the behaviour 
of as a function of Zm in the two limits, Zm — >■ Zm-i and Zm oo, we conclude that the 
solution exists if and only if the m — 1 centered scaling solution with a^j for 1 < i,j < m — 1 
given by (2.35) for /i = exists, and furthermore 

''m—1 



sign{am-i,m) = -sign ^ 



a.; 



(2.36) 



i=l 
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Finally the sign of the Hessian associated with the configuration in the /i — t- limit, after 
adding up the contribution from all critical points in the range Zj > Zi for j > i, is 



sign{arn-i,m) G{ai, ■ ■ ■ a^-i) , (2.37) 

where for 1 < z < m — 1 denote the deformed charges at = 0. Using (2.35)-(2.37) we 
get 

(m-l \ 
5^ a.,n ■ ■ ■ + J] AGfi , (2.38) 

i=l J B 

where AGb denotes the jump in G during the deformation from yU = to /i = 1 across the 
various critical points fiB where a subset B of the charges can form scaling solutions.^ Since 
the deformations involve the aim and am-3,m-i, new scaling solutions must involve either the 
mth center, or both the (m — l)'th and (m — 3)'th center. Three kinds of scaling solutions 
can be encountered during the deformation: 

1. The scaling configuration involves the charges Am- 2, and am- In this case we need 

'^m~2,m-l "I" ^ra-l,ra ^m-2,m ~ ) (2.39) 

which requires 

+ am-2,m) = . (2.40) 

2. The scaling configuration involves charges dfc, ■ ■ ■ for 2 < /c < (m — 3). In this case 
we require 

E "^. = (2-41) 

which translates to 

m—l m—1 

fiB ^aim+ E + ~ /^-b) ^aim = 0. (2.42) 

i=k i=l 

fc<i< j <m— 1 

3. The scaling configuration involves charges dfc, ■ ■ ■ «m-i for 2 < /c < (m — 3). In this case 
we require 

E = (2.43) 

k<iKj<m — l 

which translates to 

m— 1 

E &ij + (1 - I^b) E = • (2-44) 

i,i i=l 

fc<i<j<m— 1 

^Due to scale invariance a configuration where a subset of tlie centers get infinitely separated from the 
others is equivalent to a configuration where a subset of the centers come together. 
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For each of these cases the computation of AGb follows the procedure used for computing 
AFa earlier. We shall quote the final results generalizing (2.32): 

1. For fi satisfying (2.40) for m > 4 we have 

X sign {arn-l,m + «m-2,m) © ( 

(2.45) 

The case m = 4 requires special attention and will be discussed later. 

2. For fi satisfying (2.42) with k > 2 we have 



AGb = G d'l, ■ ■ ■ ^ a- x G{ai, a[ 



k-l 



i=k 



i=l 





/ 




\ 








E 


01 ij 






V 




) 


\ 



1 



5Z + 5Z 



<j <m— 1 



i=l 



The case k = 2 requires special treatment and will be discussed below. 
3. For fi satisfying (2.44) we have 



(2.46) 



m— 1 



AGb = G ■ ■ ■ d^„i, a'i, a'^ x G(dfc, d'^+i, 



sign 



X0 



E 



1 



V 



1,J 

fc<i<j<m — 1 



1 = 1 



E 



fc<i<j<m— 1 



m—1 



E 



a,; 



(2.47) 



The case where the scaling configuration involves charges a2,--- a'^ requires a special treat- 
ment. If we naively consider this as a special case of (2.46) above with = 2 or of (2.45) for 
m = 4, we would conclude that the jump vanishes since there are no scaling solution with 
two centers and hence G{a[, a'2 + ■ ■ ■ a'^) vanishes. However notice that in this case (2.41) 
with k = 2 OT (2.39) for m = 4 implies 



i=2 



(2.48) 



and hence the two centered configuration with one center with charge and the other center 
with charge + ■ ■ ■ a'^ is on a wall of threshold stability^°. As such configurations exist 



Recall that a wall of threshold stability is one on which the bound state can separate into two components 
with vanishing DSZ product. Across this wall the topology of the bound state changes but the index does 
not jump. 



- 18 - 



for vanishing FI parameters, we need to analyze the situation more carefully by working at 
fJ' = fJ'B + SfJ' where /ib is the critical value of /i at which eq.(2.41) is satisfied for k = 2. At 
this point we have 

Sfi aim for m > 4 

(d24 + A34) for m = 4 

—6uaim for m > 4 , ^ 

^ . (2.49) 

5/i (d24 + ^34) for m = 4 




In either case we can proceed to analyze the system following a similar kind of analysis used 
in computing AFa- We denote the locations of the centers as Zi, Zi = Z2 + y Wi for 2 < i < m 
with W2 = 0, Wm = 1, and look for solutions with y ~ 6fj,. The solution for Z2 and y are 
found by extremizing 

\ny (2-50) 

id 

2<i<j<m 



Wi = -^a'li 111(2:2 - zi) 



1=3 



with respect to Z2 and y, and the w^'s are given by extremizing 



Wo 



s,r 
2< s<r <m 



{2.51] 



with respect to W3, • • -Wm-i- The extremization of W2 with respect to all the Wk and the 
sign of the corresponding Hessian gives G{a2, ■ ■ ■ On the other hand the extremization 
with respect to Z2 and y gives identical conditions'^ 



Z2 - Zi 



i=3 



1 

y 



2<i<j<m 



(2.52) 



a'. 



a solution for small positive y exists for only one sign of ^ 



2<i<j<m ■' 



For small V i,j ^ -, 

— 2<i<]<m ■' 

The corresponding contribution to the sign of the Hessian can be found by taking the second 
derivative of W either with respect to y or Z2 keeping the other variable fixed, and is given 
by a multiplicative factor of 



sign 



E 



(2.53) 



11 



This can be traced to the fact that using scale invariance we can fix either Z2 or y. 
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Following a logic similar to that for AF^ and using (2.49) we find that at this critical value, 
the scaling index G jumps by 





( 


( 


\ 




agb = e 




E 




( 






\ 2<i<j<m 


I 





m—l 



Oii 



2<i<j<m — 1 



i = l 



\ 



\ 



'2' ■ ■ ■ '^m) for m > 4 



xsign(— dim) G{a'2' 
6 (-a23 («23 + «34 + "24)) X sign(a24 + "34) G(d2, ■ ■ ■ «!) for m = 4 . 



(2.54) 

This gives a recursive procedure for calculating the scaling index G(di, ■ ■ ■ dm), and therefore 
the total index F using (2.20), (2.32). The recursion is initialized by the result for three 
centers, given below. 



2.3.3 Coulomb index for 3 and 4 centers 

As a simple application of the procedure described above we shall calculate the Coulomb 
index for 3 and 4 centers. For 3 centers, collinear scaling solutions exist for 

sign(di2) = sign(d23), di3 = -di2 - d23 , (2.55) 
and the sign of the Hessian of W is (— i)9("23)+i_ xhus 

G(di, d2, ds) = e(di2d23) (-1)®("23)+1 _ (2.56) 

The total index given by (2.20), (2.32), is 

F({di, d2, da}; {gi, 52, £3}) =(-l)«(-^i2)+«(-"23) ©(^^^ g^) e(d23(ci + C2)) 

+ (_l)e{ai2)4-e(c5,3) e(di2d23) e(-(di2 + d23)(di2 + d23 + , 

(2.57) 

where the first line corresponds to Fq in (2.9) and the second line to the contribution of the 
scaling solution occurring at A = — (di2 + d23)/di3. It is straightforward, if tedious, to check 
that (2.57) agrees with the result given in [14] in a particular chamber. The result (2.57) can 
be succintly summarized by saying that F vanishes unless the sign of the 5-periodic sequence 

S123 = {ci + C2, Ci, d23, di2 + d23 + di3, di2} (2.58) 

is either constant (in which case F = 1), or flips 4 times around the sequence (in which 
case F = —1). These signs correspond to the behavior of the superpotential W at the 5 
boundaries 1/2 = oo,yi = 00, y2 = 0, ?/i = 1/2 = 0, ?/i = of the domain in which the variables 
yi, 2/2 take values. The rule (2.58) is in agreement with the existence of a gradient flow 
emanating from a critical point of W inside this domain (see Figure 1, left). 
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Figure 1: Left: The physical domain for 3-center colhnear solutions has 5 boundary components 
at which the superpotential W diverges. The sign of W on each component is that of the quantity 
indicated on the corresponding edge. Right: the physical domain for 4-center collinear scaling 
solutions also has 5 boundary components, at which the superpotential W diverges. The sign of 
W on each component is that of the linear combination of djj indicated on the corresponding edge 
or vertex. In both cases, by considering the topology of the gradient flow (indicated by the arrows 
for some suitable choice of signs on the boundary components), it is easy to convince oneself that a 
critical point exists in the physical domain if and only if the signs on the 5 boundary components 
are identical, or flip 4 times around the boundary. 

For ra = 4, we first need to compute the scaling index G{ai, ■ ■ ■ 6:4). In this case the only 
contribution to AGb in (2.38) comes from the configuration where the centers 2, 3 and 4 
come together during the deformation. Using (2.54), (2.56) we can then express (2.38) as 

G(«i, ■ ■ ■ «4) = (-l)Q("23)+e(<i34) e(ai2«23) e(-«34(«14 + «24 + «34)) 

+ (_l)e(a.3)+e(<i.4+a34) e(a23d34) e(-d23(«23 + "24 + a34)) • 

(2.59) 

The rule (2.59) can be summarized by saying that the sign of the 5-periodic sequence 

S1234 = {"12, "34, «23 + "34 + "24, "23, 012 + "23 + "13} (2.60) 

is either constant (in which case G = 1) or alternates 4 times (in which case G = —1). These 
signs correspond to the behavior of the superpotential W at the 5 boundaries Z2 = 0, = 1, 
Z2 = = 1, Z2 = Z3, Z2 = = oi the domain in which the variables Z2, z^ are valued (in 
the gauge zi = 0, 2:4 = 1). The rule (2.59) is in agreement with the existence of a gradient 
flow emanating from a critical point of W inside this domain (see Figure 1, right). 
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Using (2.20) we can now compute F{{ai, ■ ■ ■ 04}, {ci, ■ ■ ■ £4})) as a sum of four terms: 
Fo{{ai, ■ ■ ■ a^}, {ci, ■ ■ ■ £4})) given in (2.9), and the jumps across the values of /i where all 
the centers come together, where the centers 2,3,4 come together, and where the centers 
1,2,3 come together. The final result takes the form 

3 



F ({tti, ■ ■ ■ 154}, {ci, ■ ■ ■ C4}) = Y\_ 0(afc,fc+i 4)(-l)^ 



k=l 



^^e(ai3+M.4)+lC'(^«(l),...««)e (- ( Yl («12 + «23 + «34 



~(2) ~(2) ~(2) 



"2 '"3 '"4 



+ (-l)e(«24)+i^ (^1^(2)^ ~(2) ^ ~(2) ^ ~{2)^^ ^g^^ £2 + £3 + £4}) G ( 

X 0(-(a23 + «34)(a23 + "34 + «24)) 
+ (-l)e('^i3)+l^ (^1^(3) + ^(3) ^ ~(3)^ ~{3)|^ ^-^^ j ^ j^~(3)^ ~{3)^ ~ 

X e (-(ai2 + a23)(ai2 + "23 + "is)) , 



(2.61) 



where 

'^if = Ai aij for |« — j| < 2, 1 < i, J < 4, a^j'' = otherwise , 
(524'* = A2 024, tt-^'' = ctij otherwise , 

af^ = A3 ai3, tt'l'' = aij otherwise , (2.62) 

ai2 + a23 + "34 , a23 + tt34 , "12 + "23 cn\ 

Ai = —- z z — , A2 = z , A3 = z • [Z.bo) 

ai4 + ^24 + ai3 tt24 Q;i3 

This can be easily generalized to higher number of centers. 



3. Quiver invariants 

In this section we shall describe how the results of the previous sections can be used to give 
a complete prescription for computing the Poincare-Laurent polynomial of quiver moduli 
spaces. For this we need to briefly review the prescription given in [14]. 

3.1 Quiver Poincare-Laurent polynomial from Coulomb index: a review 

We shall consider a quiver with K nodes with a U (Ng) factor at the i-th node, arrows 
from the i-th node to the fc-th node representing 7^^ number of {Ni,Nk) representations of 
U {Ni) X U (Nk) and FI parameters Ci; ' ' ' Cx satisfying NiQ = 0. A negative 7^^ indicates 
7fc£ = — jik number of (iV^ , A^^) representations of U (Ni) x U {N^). Instead of considering one 
specific quiver at a time it turns out to be more convenient to consider the family of quivers 
labelled by different ranks {N^} and different values of FI parameters {Ce}- For this we assign 
to each node i a basis vector 7^ = (0, . . . , 0, 1, 0, . . . 0) - with 1 inserted at the i-th position 
- in an abstract vector space Z^', denote by F C Z-^ the collection of vectors 7 = Ylf=i ^ili 
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where Ni are non-negative integers, and by the hyperplane Yl!i=i ^id = in the space of 
real vectors ( = J2f=i Cilt ^ I^^^ • We also introduce a symplectic inner product 

K 

{-f,i)^ Y^NeK^,,, (3.1) 



£,k=l 



between the elements 7 = X]£=i of F. To any vector 7 G F and C ^ C^, we associate a 
quiver Q(7, Q with K nodes, 7^^ arrows connecting the node £ to the node fc, gauge group 
U{Ni) X U{N2) X ■ ■ ■ U{Nk), and FI parameters {Ci, ■ ■ ■ C-ft'}- If some of the A^^'s vanish we 
just drop the corresponding nodes. 

Let Q(7;C;Z/) the Poincare-Laurent polynomial 

2d 

Q{r,C;y) = J2^p{M){-yr-' (3.2) 
p=i 

where d is the complex dimension of the moduli space Ai of the quiver Q(7; C) and the 6p's 
are the topological Betti numbers of M.. The Coulomb branch formula for Q{'j',C',y), which 
we denote by Qcouiomb(7; C; v)^ takes the form: 



Qcouiomb(7;C;i/) = y2 ^ -Qcouiomh{l/m;C,y"') 

^ ^ YY) 11'''' y — '"' 

^Coul omb ({«!, ■ ■ ■ , an}, {ci, ■ ■ ■ c„}; y) 



m y^ — y~ 

m|7 

Qcoulomb(7; C; y) = / , / , ,» ,/r TNI 

^ ,ifp> |Aut({ai,--- ,aj)| 



1 - 1 

n ( E ^, ^tot(a.M; y-o) • (3.3) 



The first line is the standard relation between integer and rational BPS invariants which has 
appeared in a variety of contexts [29, 10, 30, 19]. /i(m) is the Mobius function, which is 1 
(— 1) if m is a square-free positive integer with an even (odd) number of prime factors, and 
if m is not square-free. In the second line, |Aut({ai, ■ ■ ■ «„})! is a symmetry factor given 
by Sfc! if among the set there are Si identical vectors ai, S2 identical vectors etc., 
and m\a means that m is a common divisor of (rii, ■ ■ ■ , Uk) if « = nf^n. The sums over n 
and {ai, ■ ■ ■ a„} in the second equation label all possible ways of expressing 7 as (unordered) 
sums of elements a, of F. The coefficients q are determined in terms of the FI parameters 
by Cj = A^O whenever = Ai^ii^. From the restrictions ^ ■ = 7 and NiC,i = 
it follows that Cj = 0. The Coulomb indices 5'couiomb({tti5 ■ ■ ■ , {ci, ■ ■ ■ c„}; y) can be 
computed from (2.33). The functions f2tot(a; y) are expressed in terms of the single-centered 
BPS invariants Qs through 

ntot{a-y) = ns{a;y)+ J] H{{f3i}- {m,}-y) l[ns{|3^■,y"'') ■ (3.4) 

{/3ier}.{miez} i 

Finally, the functions if({/9j}; {ki}; y) and f2s(7; y) are determined as follows. 
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1. When the number of (3i's is less that three, if {ki}; y) vanishes. 

2. For three or more number of /3j's, observe that the expression for Qcouiomb(^i ^iA; C; y) 
given in (3.3) contains a term proportional to if({/3j}; {/cj}; y) Y\i ^s(A; y^') arising from 
the choice m = 1 in the first equation in (3.3), n = 1, ai = ki(3i, nii = 1 in the second 
equation in (3.3), and mj = ki in the expression for ^totC^i ki(3i] y) in eq.(3.4). We fix 
-^({A}) {ki}jy) by demanding that the net coefficient of the product Hi ^s(A; I/''") in 
the expression for Qcou\omh{J2i ^it^i'iV) ^ Laurent polynomial in y. This of course 
leaves open the possibility of adding to if a Laurent polynomial. This is resolved by 
using the minimal modification hypothesis, which requires that H must be symmetric 
under y — )■ y^^ and vanish as ?/ — )■ oo [11]. We determine H{{(3i}] {rrii}; y) iteratively by 
beginning with the if's with three /3j's and then determining successively the ii's with 
more /3j's. 

3. H is expected to be independent of the Fl parameters and hence can be calculated for 
any value of these parameters. 

4. After determining ii({/3j}; {fcj}; y) in this way, we set f2s(7^; y) = 1 for 1 < £ < i^. For 
all other charge vectors /3, VL^{(5;y) are fixed integers, independent of y and of the FI 
parameters, which are left undetermined by the Coulomb branch analysis. Since these 
unknown constants, as well as the quivers, are labelled by the vectors a G F, there is 
one^^ unknown constant for each quiver. This can be fixed e.g. by computing the Euler 
character of the quiver moduli space for any convenient value of the FI parameters. 

As a special case of our result we can consider the case of a general Abelian quiver. This 
corresponds to 7 = Yli^nt with A^^ = or 1. As a result 7, as well as the a^'s appearing 
on the right hand side of (3.3), are primitive vectors and the a^'s are all distinct. Thus (3.3), 
(3.4) simplifies to 

n 

5'Coulomb 

({tti, ■ ■ ■ , {ci, ■ ■ ■ Cn]] y) W ^tot{oii] y) ,(3.5) 

n>l {^ier} i=l 

fitot(«;z/) = fis(«;i/)+ ^({A};y) n^s(A;z/), (3.6) 

where H{{l3i}]y) = H{{f3i}; {1,1, ■■■ l};y). The functions H{{f3i};y) are determined by 
requiring that they vanish as ?/ — )■ 0, 00, are invariant under y — y~^, and that the coefficient 
of Yli ^s{c(i] y) in the expression for (5couiomb(7; C; y) is a positive integer for each set {aj}. 

Actually the number of unknown constants is less than that of the number of quivers since ^^3(7) is non- 
trivial only if there exists a set of a^'s in F such that on — ^ and it is possible to find three dimensional 
vectors fi such that ^ oiij{ri — rj)/\ri — rj\ = 0. 
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3.2 Coulomb index for non-generic charges 

The formulae (3.3), (3.4) are completely explicit provided we have an explicit algorithm 
for computing 5'couiomb({ai, ■ ■ ■ {ci, ■ ■ ■ c„}; y). We have given such an algorithm in the 
previous sections for generic a^^s, e.g. all ajj's non-zero and no ordered subset ai, ■ ■ ■ of 
the Oj's satisfying X]i</c<f<s '^^^ ~ ^- have also assumed that the FI parameters stay 
away from the walls of marginal and threshold stability so that e.g. the quantities Yli=i ^o-{«) 
appearing in (2.10) never vanish. However we need (^couiomb for non-generic a^'s and Cj's as 
well. These come from two sources. First of all the 7^fc's of the original quiver themselves 
may be non-generic with some ^^kS vanishing or satisfying special relations. Second, even 
if the original 'jikS are generic, in the argument of (^couiomb we may have parallel ctj's. For 
these the corresponding ai/s will vanish. Also when the total dimension vector {A'^i, ■ ■ ■ A'^^-} 
is non-primitive, we shall encounter fi'couiomb in (3.3) for which the FI parameters sit on the 
threshold stability walls. In all such cases we need to evaluate S'couiomb by first deforming the 
ttjj's and/or q's to generic values and then taking the limit back to the original configuration. 
The goal of this section will be to determine a prescription for such deformations. We shall 
first describe the prescription and then justify it. 

1. To deal with the first problem we deform the •jekS to 

l£k ^ lik + ei^ik (3.7) 

where ei is a small positive number and ^^^'s are random numbers between —1 and 1 
satisfying = —^m- This will make all the generic. 

2. At this stage in any given term in the sum in (3.3), all the ctj^'s are generic except for 
subsets of Oi's which are all parallel and/or equal. In computing S'couiomb ({««}, {ci}]y) 
for a set of a^'s we consider an arbitrary ordering^^ of all the a^'s and deform them by 

ttij aij + €2 Pij, Ci ^ Ci + e2fi , (3.8) 

where €2 is a small positive number that is parametrically smaller than the previous 
paramater ei, /3ij's for i < j are randomly chosen positive numbers between and 1 with 
= and /i's are random numbers satisfying J2i fi = 0- Under such a deformation 
any subset of the {a^j's which are parallel and/or equal to each other get deformed in 
such a way that the corresponding subquiver does not contain any oriented loop. Also 
the FI parameters move away from threshold stability walls even if the undeformed 
configuration sits on such a wall. 

^■^In order to increase the efficiency of the procedure, we can arrange the nodes in some fixed order and 
then choose the ^^^'s such that > for £ < k. This will minimize the introduction of new oriented loops 
and hence scaling configurations during this deformation. For example if wc have three nodes i, j and k 
such that aij = ajk = onk = 0, then under the deformation (3.7) the subquiver containing the nodes i, j 
and k will not have any oriented loop. 

^*Hcrc we are considering the ordering as a set and not an ordering of the locations of the centers. The 
same deformation must be used for all possible arrangements of the centers. 
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3. At the end of the second step the a^j's and q's are generic and can be used to compute 
S'Couiomb- In particular the s{a) factors in (^couiomb are computed using the deformed ttj/s 
and Cj's. However in computing the ?/^«<j °'j factors in (^couiomb we use the undeformed 
ajj's, since at the end of the computation we are in any case supposed to take the aj/s 
to their undeformed values. This (^couiomb is then used to compute Qcouiomb(7; Ci v) via 
(3.3). 

In order to prove the validity of the procedure we need to argue that the deformed 
result reduces to the undeformed one in the limit when the deformations are switched off. 
First note that there is a qualitative difference between the deformations generated by ei 
and those generated by 62. For the latter the deformation of the aj/s and q's we use is 
specific to the ttj's and q's which appear in the argument of a given (^couiomb- As a result we 
need to establish that each S'couiomb returns to its undeformed value upon switching off this 
deformation. On the other hand, the deformation generated by ei can be carried out for the 
full index (5couiomb(7; Ci v) as it applies to the whole family of quivers and does not refer to 
any specific multi-centered black hole configuration. We shall indeed argue that while the 
individual (^coulomb's in the ei deformed system do not necessarily reduce to the undeformed 
result, the total index does. 

Let us begin with the €2 deformation. This contains two parts: deformation of the 
Oi/s and deformation of the q's. First consider the effect of deforming the q's. If the 
initial configuration is away from the walls of marginal and threshold stability then this 
deformation has no effect. However when 7 = ^ . is not primitve, it could happen that 
the set {ai, ■ ■ ■ a„} can be divided into two or more sets such that the sum of the a^'s in each 
set is parallel to the total charge 7. In this case the sum of the q's in each set vanishes and 
the FI parameters sit on the wall of threshold stability. The deformation of the q's given in 
(3.8) is needed to move away from this wall and make (^couiomb well defined, but the result 
does not depend on how we deform the Cj's. 

Next we turn the effect of the 62 deformation on the a^j's. We begin with the deformed 
system and take the aj/s one by one back to their values after the first deformation. During 
this process fifcouiomb can jump if two or more centers come together during the deformation. 
Using the analysis of §2.2 we know that the possible jumps in S'couiomb could arise if during 
the deformation a subset A of the centers can come in the coUinear scaling configuration by 
having tt^j = 0. Since at this stage we have already carried out the ei deformation 

making the '-fekS generic; the possible subsets where this could happen will only involve 
the centers carrying equal or parallel charges at the end of the first deformation. However 
this is ruled out by the fact that the second deformation has been choosen so that any 
subquiver, containing equal or parallel charges at the end of the first deformation, remains 
free from oriented loops. This shows that we do not encounter any collinear scaling solutions 
during the second deformation and hence there is no jump in the index Sfcouiomb during this 
deformation. 

Finally we turn to the ei deformation. To deal with this case we note that the anal- 
ysis of §2.2, showing that the refined index of a quiver changes continuously under de- 
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formations of the ai/s, breaks down on a subspace on which J2i jeA i<j '^ij ~ ^'^^ some 
subset A. Around this subspace the Coulomb index computed for a given set of charges 
will depend on the sign of the deformation parameters since a particular coUinear solu- 
tion may exist for one sign of the deformation, but as the deformation parameter ap- 
proaches zero the centers in the subset A come close together and for the opposite sign 
of the deformation parameter the solution ceases to exist. However we shall now argue 
that Qcouiomb(7; C) I/) computed from (3.4) remains independent of the sign of the defor- 
mation. For this suppose that we have three centers carrying charges ai, 02 and such 
that ai2 + + ctis = 0, ai2,«23 > 0. In that case for the permutation (123) we can 
have a coUinear scaling solution. Now in the deformed system whether this permuta- 
tion contributes or not depends on whether au, 023 and 031 form an oriented triangle 
or not, and this in turn will depend on the details of the deformation. The difference in 
fi'Couiomb({tti5 C(2, ^s}; {^1, C2, C3}; y) that we shall get between these two cases is, up to a sign, 
given by {y — y~^)~'^ since y^i<j°''^^ = 1. In the expression for Qcouiomb(«i + 02 + «3; Ci v) 
there will be a term proportional to 5'couiomb({ai, "2, as}; {ci, C2, C3}; ?/)fis(ai)fis(a2)^^s(a3) 
and the ambiguity in g'couiomb described above will lead to an ambiguity proportional to 
(y - y"^)"2fis(ai)^s(a2)fis(a3) in the expression for Qcouiomb(ai + "2 + a3;C;i/)- Now 
the expression for Qcouiomb(ai + ^2 + ct-s'Xiy) "will also contain a term proportional to 
H{{ai, a2, a-s}] y)Qs{c(i)^s{c(2)^s{c(3) , and the function H{{ai, 02, 03}; y) is determined by 
requiring that H vanishes as y — t- 0, 00 and that the net coefficient of f2s(cn)fis(a2)^s(a3) 
is a polynomial in y, y~^. Since {y — y"^)'^ — )■ as ?/ — )■ 0, 00 we see that the ambiguity in 
5'Couiomb introduced above is absorbed completely into the function H and does not lead to 
any ambiguity in the expression for Qcouiomb(ai5 ^25 ^3; C; ?/)• This argument can be easily 
generalized to argue that all the ambiguities in (^couiomb can be absorbed into the functions 
H and the expression for Qcou\omb{l', Ci u) is independent of the choice of the deformation of 
the 7ffc's for general vector 7. 

Using the analysis above we can also answer the question: how small should ei and €2 
be for the procedure described above to hold? The general rule is that if a'^j denotes the 
deformed aij, then during the deformation we should not have ^i^j-i j^^A^'ij ~ 
ordered subset A of the centers. This means that we should not encounter any coUinear 
scaling configurations during the deformation except possibly at the beginning. 

4. Coulomb/Higgs equivalence for quivers without loops 

For quivers without oriented loops, quiver invariants can be computed using the Harder- 
Narasimhan recursion, or equivalently using Reineke's solution to this recursion [15]. In this 
section we shall show that for such quivers the Coulomb branch formula (3.3) agrees with 
Reineke's formula, both for Abelian and non-Abelian quivers. 

^"''This would happen if the mitial configuration had J2i<j-i j^A'^ij ~ ^ for some ordered subset A of the 
centers. 
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4.1 Reineke's formula for quivers without loops 



As reviewed in [14], the Poincare-Laurent polynomial of a quiver Q without oriented loops 
can be computed using the Harder- Narasimhan recursion. Henceforth we shall denote the 
Poincare-Laurent polynomial computed by this method by Qmggs{^£Ni'yfX;y), to distin- 
guish it from the Coulomb branch formula (3.3). The expression for Qmggs takes the form 
[12], [16, Theorem 6.8] 

QHiggs(7; C; = 2^ ^ ^^_^-,n ^Higgs(7M; C; y ) 

m|7 

QHiggs XI ^' ^ = X] 



l(y-\lyY-^ (4.1) 

Y}- , A/(*)=JV;J\7{*) IliVfor i = 



X n GHiggs({M«, ■ ■ ■ M«}; {71, • ■ • IkY {Ci, " " " U}\ y) 

i=l 



The sum over {M^*^} runs over all orderec? partitions of iV = (A^^i, ■ ■ ■ Nk) into parallel vectors 
. . . and Gmggs{.{M^\ ■ ■ ■ M^^}; {71, ■ ■ ■ 7^}; {Ci, ■ ■ ■ Ck}] y) is the 'stack invariant' 
associated to the quiver with dimension vector M*^*^ If M^*) is primitive then the quiver 
moduli space is smooth and Gniggs is just its Poincare-Laurent polynomial, however in general 
it is a rational function of not necessarily invariant under y ^ \/y. In all cases however, 
we assume that it is given by Reineke's formula [15],^^ 

GHiggs({Mi, ■ ■ ■ Mi}; {ai, ■ ■ ■ a^; {ci, ■ ■ ■ cz.}; y) = (-y)"^'^^ a/,m, max(a.„o)-i+E.M. 

X (y2_l)l-E,M, (-l)^-ly2E„<5E,,maxK,,0)^,^JV; -QQ^a^^jl^-l^ (4.2) 

partitions a,i 



where the sum over partitions in (4.2) runs over all orfiered partitions of the vector (Mi, ■ ■ ■ Ml) 
into non-zero vectors {{N^, ■ ■ ■ N^), a = 1, . . . , s} for s = 1, . . . , Mj, satisfying A''"' > 0, 
J2a = (assuming that ^ • MjCj = 0) 

b L 

E > (4-3) 

a=l i=l 

for all h between 1 and s — 1. [A^, y]! denotes the g-deformed factorial, 

[N, y]\ ^ [1, 2/][2, y]... [AT, y] , [A^, y] ^ . (4.4) 

- 1 



'It is related to the quantity 1(7, w) defined in [14], Eq. (2.39) via GHiggs(7; C; 2/) = iX/v ^ v) ^v)- 
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It is worth noting that (4.1) can be inverted to express the stack invariants in terms of the 
rational Poincare-Laurent polynomials [16, Theorem 6.8]: 

CffiggsliiVi, ■ ■ ■ Nk}; {7i, ■ ■ ■ Ik}; {Ci, ■ ■ ■ Ck}; y) 

= Y1 Yl ^y7^73^;3ninn^HisssK;C;i/), (4.5) 

k {a,} '^^ ^ ' i=l 

Eti "i=E^ Nae 

"illEiv^^^ for i=l,...,k 

where again the sum over runs over all ordered partition of Yle Ni'^i into parallel vectors 
«!, Q!2, ■ ■ ■ ttfc- We shall refer to the formula (4.1) (or equivalently (4.5)) and (4.2) as the Higgs 
branch formula for the Poincare-Laurent polynomial Q(7; C; ?/)• Our goal will be to show the 
equality of this Higgs branch result with the Coulomb branch formula (3.3). 

If Ni'-)i is a primitive vector, i.e. the iV^'s have no common factor other than unity, 
then (4.1) reduces to a simple form 

Qffiggs Nau C; y^ = GmsUi^u ■ ■ ■ Nk}; {71, ■ ■ ■ 7x}; {Ci, ■ ■ ■ Ck}; y) ■ (4.6) 

In this case Gniggs is a symmetric Laurent polynomial since Qmggs is. 
4.2 Abelian quivers without loops 

For quivers with dimension vector Ni = 1, the stack invariant Gniggs appearing in (4.6) is 
given by (4.2) where the integers A'^'^ can only be equal to or 1. We shall use a special 
symbol (^Higgs for labelling the corresponding Gniggs^ 

^Higgs({7i, ■ ■ ■ Ik]] {Ci, " " " Ck]] y) = GHiggs({l, ■ ■ ■ 1}; {71, ■ ■ ■ Ik]] {Ci, ■ ■ ■ Ck]] y) ■ (4.7) 

Assuming for the moment that all 7jj's are non- vanishing, we shall choose a strict ordering 
convention for the 7j's such that 7^ > iff i < j. Thus (4.6), (4.2) can be expressed as 

QHiggs(7i H ^lK]C]y)= 5'Higgs({7i, ■ ■ ■ Ik]] {Ci, " " " Ck]] y) = 

(^y _ y^iy-K (_ 1)^-1^2 E„<6E,<i7,«^fiv}-Ei<i<,<K 7., . (4.8) 

partitions 

where the sum runs over all ordered partitions of 7 = 71 + - ■ ■+'Jk into vectors (3^°"^ = Niji 
with A^" = 0, 1, a = 1, ■ ■ ■ s, satisfying 

b n 

E ^^C. > . (4.9) 

a=l i=l 

Since (4.8) is continuous at 7jj = 0, the result when some of the ji/s vanish can be obtained 
as a limit of (4.8). Furthermore, throughout this subsection we shall work with generic FI 
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parameters for which the left hand side of (4.9) never vanishes. In that case (4.8) is also 
invariant under small deformations of the FI parameters. 

On the Coulomb branch side, the Poincare-Laurent polynomial Qcouiomb(7; C; y) coin- 
cides with its rational counterpart Qcouiomb(7; Ci 2/)- Since the quiver has no oriented loop 
there are no contributions from scaling solutions, therefore we can set Qtot{c(,y) = ^s('^;2/) 
in (3.5). We can further set ^3(0^; u) to except when a is equal to one of the basis vectors 
7^ in which case its value is 1. The Coulomb branch formula (3.5) thus reduces to 

Qcoulomb(7l H ^lK;C,y) = fl'Coulomb({7l, " " " Ik}] {Cl, " " " Ck}] V) ■ (4.10) 

It follows from the discussion in §2.2 that the result when some of the 7ij's vanish can be 
obtained as a limit of (4.10) for generic non-vanishing 7ij's. Comparing (4.10) and (4.8) we 
see that the proof of equivalence of Qcouiomb and Qniggs for Abelian quivers reduces to the 
proof of equivalence of (^couiomb and gmg^s- Furthermore since for both (^couiomb and gmggs the 
result when some of the 7ij's vanish can be obtained as limits of the result for generic 7ij's, 
it is enough to prove the equivalence for generic non- vanishing 7ij's. 

This equivalence can now be proved by following exactly the same analysis as §3.3 of 
ref.[10]. We summarize the main points, refering the reader to [10] for more details. First, 
to each permutation a, we associate a family of ordered partitions of 71 -|- ■ ■ ■ 7^ into vectors 
as follows: 

• Break the sequence {o"(l), cr(2), ■ ■ ■ cin)} into subsequences {cr(l), cr(2), • • ■ cr(ii)}, {a{ii+ 
1), ■ ■ ■ (7(22)}, ■ ■ ■ , {o"(«a-i + l), • • • , cr(^a)}, ■ • ■ , Q = Iq < ii < • • • < is = K ^ such that each 
subsequence represents an increasing subsequence: a{ia~i + l) < a{ia-i+2) < ■ ■ ■ < a{ia) 
for each a. s gives the number of such increasing subsequences. 

• For each such breakup of {cr(l), ■ ■ ■ o-{K)} into increasing subsequences, we can associate 
a partition of 71 -|- ■ ■ ■ 7^- into vectors f^^""^ as follows: 

[3^'^'-''^= g 7.«- (4.11) 

This generates an ordered partition of 7 but does not necessarily satisfy (4.9). It is easy 
to show that with this choice the power of y in y^i<j ^<^(i)'^u) in (1.9) matches the power 
of y inside the sum in (4.8) 

For illustration we can consider the permutation (t(1234) = (3412). The increasing 
subsequences are {{34}, {12}}, {{3}, {4}, {12}}, {{34}, {1}, {2}}, and {{3}, {4}, {1}, {2}}. 
Associated partitions are {71+72,73+74}, {71+72,74,73}, {72,71^73+74} and {72,71^74,73}, 
respectively. For each of these partitions we have 

2 ^ ^ ^jiN^Nll - -fij = la{i)a{j) = 712 + 734 " 7l3 - 7l4 " 723 " 724 • (4. 12) 

a<b j<i l<i<i<4 1<«<J<4 
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The above analysis shows how to each permutation a we can associate a set of ordered 
partitions of 71, ■ ■ ■ into vectors (3^^\ ■ ■ ■ f3^^\ The converse is also true - all ordered parti- 
tions of the vectors (71, . . . , 7i^), before imposing (4.9), are in one to one correspondence with 
the set of all increasing subsequences of all the permutations of (12 . . .n) via the rule (4.11). 
For a given permutation we shall call an increasing subsequence maximal if it is not possible 
to build bigger increasing subsequences involving the elements of the subsequence. In the 
example described above the maximal increasing subsequences are {{34}, {12}}. The com- 
plete contribution to the coefficient of |/^^<i '''■'(O'^Cj) in the Reineke formula can be generated 
by beginning with the maximal increasing subsequences associated with the permutation a 
and combining them with the contribution from other increasing subsequences associated 
with the same permutation, but we must be careful to pick only those partitions which 
satisfy (4.9). As shown in [10], a given permutation contributes if and only if its maximal in- 
creasing subsequences generate a partition via (4.11) satisfying (4.9), and none of the other 
(non-maximal) increasing subsequences generate an allowed partition. The sign of the 
contribution is given by (— l)'*"^ where s is the number of maximal increasing subsequences. 

It is now easy to see that the condition that the maximal increasing subsequences gen- 
erate a partition satisfying (4.9) translates to 

K 

Ca{€] > for 7a(fc),<7(fc+i) < 0, (4.13) 

i=k+l 

and the condition that none of the non-maximal increasing subsequences satisfy (4.9) trans- 
lates to 

K 

c<^(i) < for 7a(fc),<7(fc+i) > 0. (4.14) 

i=k+l 

These precisely agree with (2.7). Finally the number of maximal increasing subsequences is 
one more than the number of negative 'ya-(k),a(k+i)- Thus we have 

K-l 

i-iy-' = Yl sign(7<,(fc),<,(,.+i)) , (4.15) 

k=l 

in agreement with (2.8). This proves that the Reinecke formula (4.8) agrees with the Coulomb 
index (1.4), and hence also the Coulomb branch formula (4.10) for the Poincare-Laurent 



^^In particular, if some non- maximal increasing subsequence satisfies (4.9), then the net contribution from 
the partitions associated with the maximal and non-maximal increasing subsequences of a given permutation 
cancel pairwise. 
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polynomials^, 

fl'HiggsllTi, ■ ■ ■ Ik}] {Ci, ■ ■ ■ Ck}; y) = fi-couiombllTi, ■ ■ ■ Ik}; {Ci, ■ ■ ■ Ck}; y) , 
QHiggs(7i H ^iK]C\y) = Qcouiomb(7i H ^iK;C,y) ■ 

(4.16) 

For charge configurations relevant for wall-crossing, where all the 7^'s lie in a plane and the 
(iS are determined in terms of jij, this equivalence was proved by induction in [28]. Here 
we have shown that it holds for any Abelian quiver without oriented loop. 

4.3 Non-Abelian quivers with primitive dimension vector and without loops 

For non-Abelian quivers without loops we still can set Qtot{ci',y) = ^s{ci',y)- We also have 
fls{ci',y) = unless a = ■ji, and Qsili'jy) = 1- Thus r2s(Q!i/mj; i/™') is non-vanishing (and 
equal to 1) only when ai/rrii is a basis vector. Furthermore for primitive charge vector, 
<5couiomb(7; C; y) is equal to Qcouiomb(7; C; y)- in that case we can express (3.3) as 

^7couiomb({(^7,)'^''};{(^0)'^"};2/) 

^{nn^(^)'"}. (-^) 

The argument of fi'couiomb given above corresponds to a choice of {ai, ■ ■ ■ a„} where we have 
fcj^'' copies of the vector i'jj for ^ = 1, 2, 3, ■ ■ ■ , j = 1, ■ ■ ■ i^, and the FI parameter associated 
with each such copy is iQ. It is understood that if some of the Oj/s appearing in the 
argument of f^couiomb vanish then we must deform them in a way so that the deformed quiver 
does not contain any oriented loop, compute (^couiomb and then take the deformations back 
to zero. The Yli Ylt ^i^^' factor in the denominator represents |Aut({ai, ^2, ■ ■ ■ , an})|- 

On the other hand, it was shown in [10, 21] that the stack invariants (4.2) satisfy the 

^^This relation holds away from walls of marginal and threshold stability. On such a wall, the Higgs 
branch index on the l.h.s. given by (4.8), (4.9), is still well defined but the derivation of the Coulomb index 
given in (2.7), (2.8) clearly breaks down. Furthermore even the Higgs branch index given by (4.8), (4.9) 
does not have the form of a symmetric Laurent polynomial, and hence cannot be interpreted as Tr'{—y)^^'^ 
associated with a quantum system. On a threshold stability wall we shall define gcouiomb by deforming the 
FI parameters away from the wall since there is no independent definition of 5couiomb in this case, and the 
result is independent of the deformation. However for gniggs we shall continue to use the definition given in 
(4.8), (4.9) for reasons which will become clear in §4.4. The price we pay is that gcouiomb and ffHiggs are no 
longer equal on the thereshold stability wall. 
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Abelianization formula^^ 

Gm,U{NjhHh{Qhy)= E G'Higgs({iv,};a};{0};y)n 



h\ \^{y^ - y-^) 



where, on the r.h.s, GHiggs({^j}; {7j}; {Cili v) is the stack invariant of a quiver Q defined as 
follows: 

• The vertices of Q are obtained by replacing the vertex i of Q by a collection ^ of 
vertices with k = 1, . . .ki for any ^ in the partition Ni = "^^iki, and keeping all the 
other vertices j ^ i of Q; 

• Each arrow j ^ k in Q with j,k ^ i induces an arrow j — )► A; in Q; 

• The arrows i j in Q (resp. j — > i) with j i induce i arrows ig^k — > j (resp. j — )■ i^^fc) 
in Q for each i, k; 

• The dimension at each node ii^k in Q is equal to one, hence justifying the name; the 
dimensions at the other nodes j i are the same as the dimensions Nj in the original 
quiver Q; 

• The FI parameter at each node ie^k in Q is equal to £Q] the FI parameters at the other 
nodes j ^ i are the same as the FI parameters (j in the original quiver Q. 

Thus for example on the right hand side of (4.18) the set {A'^} will include all the iV^'s 
except Ni and J2e^f^ number of I's, the set {^j} will include all the 7/s except 7^, ki copy 
of 7i, /c2 copy of 2'ji, etc., and the set {Q} will include all the C/s except Q, ki copies of 
Ci, ^2 copies of 2(i etc. By successive application of this formula we can express the stack 
invariant of any non-Abelian quiver in terms of that of a family of Abelian quivers: 

K /,_,-!■ 

GHiggs({A^.}; {7^}; {0}; y) = E^Higgs({(^7,)'^''}; {(^0)''''}; y)UU ^ ^ 



#! V{y'-y~') 



(4.19) 

Using (4.6) we can replace the left hand side by Qniggs {^f=iNili'-,C'iy^ ■ The resulting 
equation is identical to (4.17) with Qcouiomb replaced by Qniggs and (^couiomb replaced by 
S'Higgs- Using the equivalence of fi'couiomb and gm^gs proven in (4.16), we get 

Qffiggs N^t] C; y j = Qcouiomb N^f, C; yj ■ (4.20) 



^^In [10], appendix D we proved this formula for a special choice of FI parameters relevant for wall 
crossing. However this assumption was inessential and the result easily generalizes to the case of general FI 
parameters [21]. The formula (4.18) also holds in the case of quivers with loops, but wc shall not make use 
of this fact here. 
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Thus the Higgs and Coulomb branch formulae are equivalent for non-Abelian quivers as well, 
so long as the dimension vector {Ni} is primitive. 



4.4 Non-primitive dimension vector 

We now turn to the case of a non-Abelian quiver with non-primitive dimension vector. In 
this case the Coulomb branch formula is still given by (3.3) although Qcouiomb and (Qcouiomb 
are no longer identical. On the other hand the Higgs branch formula is given by (4.1) (or 
equivalently (4.5)) and (4.2). Our goal in this section will be to prove the equality of Qcouiomb 
and Quiggs- 

When the dimension vector {Ni}i=i_,,K = N is not primitive, the Abelianization formula 
(4.19) is still known to hold, and can be used to express Gmggs in terms of 5'Higgs({(-^7j) ^ }; 
{(^0)*^^ };?/) for integers kj^'^ satisfying Yli^^f' ~ ^j- However some of the Abelian stack 
invariants 5'Higgs({7j}; {0};?/) appearing on the r.h.s. have stability conditions {(} lying on 
walls of threshold stability. This happens when the set {7j} can be divided into several 
subsets Ai, ^2, ■ ■ ■ such that J2ieAs ^« II ^^^^ ^^^^ ^- ^ resuh 5'Higgs({7i}; {Ci}] v) 
can no longer be equated to 5'couiomb({7i}; {0}; u) (see footnote 18). Let {Q}' be a sufficiently 
generic perturbation of the FI parameters away from the wall and sufficiently small so that 
no other walls of marginal stability are crossed. In this case 5'Higgs({7i}; {0}'; v) can again be 
equated to 5'couiomb({7«}j {0}'iZ/) and furthermore will be independent of the deformation. 
Now the HN recursion method [31, 15] relates 5'Higgs({7i}^ {C«}') v) and 5'Higgs({7i}; {Oji v) as 
follows: 

5'Higgs(«l, ■ ■ ■an;Cl, ■ ■■Cn]y) 

k 

= Y.^y-'-yf~' E n^Higgs({a.; ^^A^■{d^ zG^};y) (4-21) 

fe>l {As} s=l 

where the sum over {A^} runs over all unordered partitions of the integers 1, 2, • • ■ n into sets 
Ai, ■ ■ ■ Ak subject to the conditions indicated above. The symbol {c^} denotes a generic 
deformation in which the FI parameters associated with each of the n nodes are deformed 
independently so that the configuration moves away from the wall of threshold stability, 
subject to the constraint that the sum of the FI parameters carried by all the nodes must 
vanish. It was shown in Ref. [22] that the combinatorics of the summation in (4.21) can be 
summarized by an equality of generating functions, which in our notation reads 



1 + iy-' - yr' E ^Higgs ({(^7,)'^^''}; m,)'^''}; y) U 



it 



ii 



:exp [y-'-y)-' ^ ^^mggs ({(^7,)'^"}; {(^0)'^^''}'; Z/) H 



(4.22) 
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where tj^ are formal parameters and the sum over {k^^^} on either side runs over all integers 
kf"^ for which A;f^^ is of the form cNi for some fixed vector Ni and arbitrary constant c. 
Note that for the deformed stabihty conditions (^Higgs can be equated to S'couiomb, which 
is by definition equal to the value of (^couiomb for the undeformed stability condition ( (see 
footnote 18). As a result, we get for (^Higgs either from (4.21) directly, or after equating 
coefficients of Ylii{ti,e)'''- on either side. 



l-k 



il 



k\ 



X 



n |^?couion.b ({(^7,)'^'^'}; {(^0)'^'^'}; 2/) n ij 



(4.23) 



where the sum over {kf'^'^} runs over all ordered partitions of {kf^}. Using this we can 
express the Abelianization formula (4.19) as 



G'Higgs({A^.};{7^};{C};y) 



E (n*!"OE 



^ 3 ' 



k\ 



E 



{fc(^'''',l<s<fe} 



n 

s=l 



^7couion.b({(^7,)'^"^'};{(^0)'^'"};y)n 



ii '^i 



K 



nn 



(4.24) 



We can now remove the sum over {kf''} by relaxing the constraint kf'^^ = k\''\ replacing 
kf^ by fcj-^'*'' everywhere, and imposing the constraints ^ ikf'^'^ 
Ni. This gives 



A^i and Y.i^kf''^ « 



Gmggs({Ar.};{7.};{0};y) 

k 



E 



A;! 



E 



{kV-''\l<3<k} 



n 

s=l 



^?couiomb {(^7,)'^ };m,P hy 



X 



y-y 



-1 



(4.25) 



Following the discussion at the beginning of §4.3 but without imposing the primitivity 
condition on the dimension vector 7 = A£7£, we see that the Coulomb branch formula 
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(3.3) takes the form 

<ycouiomb(7; C;y> = -<ycouiomb(7/"^; C; v > 

Z — / jji ym _ y-m 

m|7 

K 1 / _ -1 N 

gcouiomb(7; C; y) = E ^?couiomb({(^7,)'=^''}; {(^0)'^''}; ^) H H 7w7 pL -y-e ) 

^ 3 ' 

(4.26) 

Let us now define Gcouiomb by a formula analogous to (4.5) with Qniggs replaced by Qcouiomb 
on the right hand side: 

Gcouiomb({iVi, ■ ■ ■ Nk}] {7i, ■ --Ik}] {Ci, ■ ■ ■Cx};^) 

= Y. kUv - n Ocouiomb(«^; C; y) ■ (4.27) 

Proving the equivalence of Qniggs and Qcouiomb is then equivalent to proving the equivalence 
of Gniggs and Gcouiomb- Now using (4.26) we get 

Gcouiomb({iVi, ■ ■ ■ Nk}; {71, ■ ■ ■ Ik}] {Ci, ■ ■ ■ Ck}] v) 

= E E Ffi^^^n E ^co.io.b({(^7.)^'''};{(^0)^''"};y) 



X 



1 f y-y ' ^ ' 



nn 



(4.28) 

where the sum over {k^^'^^} runs over all or(ierec? partitions of {kf^}- We can now remove the 



sum over a^'s by relaxing the constraint that ^^^k'^j'^'^lj = c^s, but imposing the constraints 
Y.^^Jkf'''^ = Ni and E^^^f''^ « iV^. This gives 

G'couiomb({iVi, ■ ■ ■ Nk}] {71, ■ ■ ■ Ik}] {Ci, ■ ■ ■ Ck}] y) 

k 

=T.l_(y''-yy~' E E ^coniomb({(^7.)'^'^'};{(^o)'^'^'};y) 

^ {fc(^'^\l<s<fc} - 

L J 1 _ _ J 



X 



nn 



(4.29) 
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Since the right hand side matches the right hand side of (4.25), it follows that Gniggs and 
Ccouiomb are identical, and therefore the Higgs and Coulomb branch computations Qmggs 
and Qcouiomb are equivalent, even for non-primitive dimension vector. 

Given the equivalence of Qmggs (defined by the Harder-Narasimhan recursion (4.1), (4.2)) 
and Qcouiomb (defined by the Abelianization formula (4.26)), and given the equivalence of 
the Abelian indices (^couiomb and (^Higgs for deformed FI parameters, we conclude that the 
Poincare-Laurent polynomial Qniggs of a quiver without loop satisfies, for arbitrary dimension 
vector, the Abelianization formula: 

Qffiggs(7; C,y)=2^ — ^;;73^QHiggs(7M; C; y ) 

m|7 

Qmggs(7; C,y)= E ^7mggs({(^7,)'=''}; im'^'}'; ^) H H ^ [Ju^^ 



(4.30) 



This can be regarded as the main result of this subsection. 



A. The Mathematica package "CoulombHiggs" 

To facilitate further investigation, we provide a Mathematica package allowing to com- 
pute the Coulomb index (^couiomb for multi-centered black holes and the Poincare-Laurent 
polynomial Qcouiomb(7; Ci I/) and QHiggs(7; Ci I/) of quiver moduli spaces using the Coulomb 
branch and Higgs branch formulae. We also provide three example files where this package is 
used to evaluate quiver invariants for the Kronecker quiver (Kronecker .nb), for non- Abelian 
3-node quivers (Threenode .nb) and for several 4 and 5-node Abelian quivers considered in 
[14] (Multinode .nb). The validity of the algorithm for the F and G indices is tested in a 
fourth file CoulomblndexCheck.nb. All these files are included in the "source" of this paper 
available from arXiv and can be obtained from the second name authors' webpage. 

Assuming that the file CoulombHiggs .m is present in the user's MATHEMATICA Appli- 
cation directory, the package is loaded by entering 

iii[i]:= <<CoulombHiggs ' 

out[i]:= CoulombHiggs v 1.0 - A package for evaluating quiver 

invariants using the Coulomb and Higgs branch formulae. 

If the file CoulombHiggs .m has not yet been copied in the user's Mathematica Appli- 
cation directory but is in the same directory as the notebook, evaluate instead 

in[i]:= SetDirectory [NotebookDirectory [] ] ; <<CoulombHiggs ' 
out[i]:= CoulombHiggs v 1.0 - A package for evaluating quiver 

invariants using the Coulomb and Higgs branch formulae. 
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The first main routine is CoulombBranchFormula, whose basic usage is illustrated below: 



in[i]:= Simplify [CoulombBranchFormula [4{{0, 1, 0, 1}, {1, 

-1, 0}}, {1/2, 1/6, -2/3}, {1, 1, 1}, y]] 

Out[l]:= 2+ ^+?/2 + 0mS({l,l,l}) 

The first argument corresponds to the matrix of DSZ products aij (an antisymmetric 
matrix of integers), the second to the FI parameters Q (a vector of rational numbers), the 
third to the dimension vector Ni (a vector of integers), and the last to the angular momentum 
fugacity (a formal variable). In this example, the routine expresses the Poincare-Laurent 
polynomial of the moduli space of a three-node Abelian cyclic quiver with 4 arrows between 
each subsequent node, in terms of the single centered degeneracy fig (71 + 72 + 73). 

The second main routine is HiggsBranchFormula, which computes the Poincare-Laurent 
polynomial using the Higgs branch formula (4.1). The arguments are the same as for 
CoulombBranchFormula, but it only applies to quivers without oriented loop (although no 
test is made to ensure that this is the case): 

in[i]:= Simplify [HiggsBranchFormula [{{0, 3}, {-3, 0}}, {1/2,-1/2}, 
{2, 2}, y]] 

Out[l]:= F 

The above command computes the Poincare-Laurent polynomial for the Kronecker 
quiver with 3 arrows, FI parameters (1/2,-1/2), dimension vector (2,2). The package 
allows for much more, as documented below. Inline documentation can be obtained by 
typing e.g. 

in[i]:= ?CoulombBranchFormula 

Out[l] : = 



A.l Symbols 

• y: Angular momentum fugacity; 

• Om[ci3arge yector_,y_ ]: denotes the refined index ^(7,?/); 

• Omb[ci2arge yector_,7_ ]: denotes the rational refined index 0(7,?/); 

• OmS[ci2arge yector_,j_ ]: denotes the single centered degeneracy fis(7; u)- The dependence 

on y can be omitted, but only at the end of the computation; 

• OmT[ciiarge vector_,7_ ]: denotes the (unevaluated) function f2tot(7,l/) defined in (3.4); 

• Coulombg[]ist of charge vectoTs_,y_ ] : : denotes the (unevaluated) Coulomb index (/couiomb ({««}, {q}, y) 

leaving the FI parameters unspecified; 
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• HiggsGCciiarge yector_,y_ ]: denotes the (unevaluated) stack invariant Gmggs{l,y) defined 

in (4.2); 

• CoulombH[]ist of charge vectors_,multiphcity vector_,y_'\: denotes the (unevaluated) factor 

H{{ai}, {ni},y) appearing in the formula (3.4) for Qtot{^niai,y) in terms of f2s(ai,y)- 

• QFact[ii_,7_]: represents the (non-evaluated) g-deformed factorial [n,y]\ 
A. 2 Environment variables 

• SQuiverPerturbl: Sets the size of the perturbation ei = l/$QuiverPerturb of the DSZ 

products in (3.7), set to 1000 by default. 

• $QuiverPerturb2: Sets the size of the perturbation £2 = l/$DSZPerturb of the DSZ 

products in (3.8), set to 100000 by default. 

• $QuiverNoLoop: If set to True, the quiver will be assumed to have no oriented loop, hence 

all H factors and all fis('^) "will be set to zero (unless a is a basis vector). Set to False 
by default. 

• SQuiverTestLoop: If set to True, all H factors and f2s(ct) corresponding to subquivers 

without loops will be set to zero (unless a is a basis vector). Set to True by default. 
Determining whether a subquiver has loops is time-consuming, so for large quivers it 
may be advisable to disable this feature. Note that SQuiverNoLoop takes precedence 
over this variable. 

• $QuiverMultiplier: Overall scaling factor of the DSZ matrix in any evaluation of Coulombg 

or HiggsG. Set to 1 by default, could be a formal variable. 

• $QuiverVerbose: If set to False, all consistency tests on data and corresponding error 

messages will be skipped. Set to True by default. 

• SQuiverDisplayCoulombH: If set to True, the routine CoulombBranchFormula will return 

a list {Q, R} where Q is the Poincare-Laurent polynomial and i? is a list of replacement 
rules for the CoulombH factors. Set to False by default. 

• $QuiverPrecision: Sets the numerical precision with which all consistency tests are 

carried out. This is set to by default since all data are assumed to be rational numbers. 
This can be set to a small real number when using real data, however the user is warned 
that rounding errors tend to grow quickly. 
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A. 3 Coulomb index 

• CoulombF[Mat_,Cycc_] : returns the index of collinear solutions F{{ai, ■ ■ ■ {ci, ■ ■ ■ c^}) 

with DSZ products aij = Mat[[i, j]], Fl terms = Cvec[[i]] and trivial ordering. 

• CoulombG[Mat_]: returns the index of scaling collinear solutions G{{ai, ■ ■ ■ an}) with 

DSZ products aij = Mat[[2,j]] and trivial ordering. The total angular momentum 
Mat[[2, j]] must vanish; 

• CoulombIndex[Mafc_,PMafc_, Cvcc_,y_ ] : evaluates the Coulomb index (^couiombllcti, ■ ■ ■ ctn}; 

{ci, ■ ■ ■Cn};y) with DSZ products a^j = Mat[[i,j]], perturbed to PMat[[i,j]] so as to lift 
accidental degeneracies, possibly rescaled by an overall factor of SQuiverMultiplier, 
FI terms Cj = Cvec[[i]], angular momentum fugacity y; 

• CoulombFNum[Mat_] : computes numerically the index F{{ai, . . . an}, {ci, . . . c„}) with DSZ 

matrix ctij = Mat[[i, j]] and FI parameters Cj = Cvec [[?]]. For testing purposes only, works 
for up to 5 centers. 

• CoulombGNuni[Mat_]: computes numerically the scaling index G{ai, . . . an) with DSZ ma- 

trix aij = Mat[[i, j]]. For testing purposes only, works for up to 6 centers. 

• CoulombIndexNum[Mat_,PMat_, Cvec_,if_,j_ ] : returns the Coulomb index gcouiomh{{c(i, ■ ■ ■ 

an}] {ci, ■ ■ ■Cn};y) with DSZ products aij = Mat[[z,j]], possibly rescaled by an overall 
factor of SQuiverMultiplier, FI terms q = Cvec[[z]], angular momentum fugacity y, 
by searching collinear solutions numerically; For testing purposes only, works for up to 
5 centers. 

A. 4 Coulomb branch formula 

• CoulombBranchForniula[Mat_, Cvec_,iVyec_,.7_ ]: computes is the Poincare-Laurent poly- 

nomial of a quiver with DSZ products aij = Kat[[i,j]], dimension vector iV, = Nvec[[i]], 
FI parameters Q = Cvec [[«]], in terms of single-centered invariants fls- This stan- 
dalone routine first constructs the Poincare-Laurent polynomial using (3.3), evaluates 
the Coulomb indices (^couiomb, and determines the H factors recursively using the min- 
imal modification hypothesis. If SQuiverDisplayCoulombH is set to True, the routine 
returns a list {C|,R}, where Q is the Poincare polynomial and R is a list of replacement 
rules for the CoulombH factors. For quivers without loops, the process can be sped up 
greatly by setting SQuiverNoLoop to True. For complicated quivers it is advisable to im- 
plement the Coulomb branch formula step by step, using the more elementary routines 
described below. 

• CoulombBranchFormulaFromH[Mat_, Cyec_,iVvec_,P_,7_ ]: returns the Poincare-Laurent 

polynomial of a quiver with DSZ products aij = Mat[[2,j]], dimension vector Ni = 
Nvec[[z]], FI parameters Q = Cvec[[z]], using the rule R to replace all CouiombH factors. 
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• QuiverPoincarePolynomial[iVyec_,y_]: constructs the Poincare-Laurent polynomial of a 

quiver according to (3.3). Coincides with QuiverPoincarePolynomialRat for primitive 
dimension vector; 

• QuiverPoincarePolynoinialRat[iVvec_,y_] : constructs the rational Poincare-Laurent poly- 

nomial (5couiomb(7;C;i/) according to (3.3); 

• QuiverPoincarePolynomialExpand[Mat_,PMat_,Cyec_, Nvec_, Q_]: evaluates the Cou- 

lomb indices (^couiomb and total single center degeneracies Qtot{cii,y) appearing in the 
Poincare-Laurent polynomial Q of a quiver with DSZ products aij = Mat[[i, j]], per- 
turbed to PMat[[z, j]], dimension vector Ni = Nvec[[i]], FI parameters Q = Cvec[[i]], 
using (1.9) and (3.4); 

• CoulombHSubQuivers[Mat_,PMat_,iVyec_,y_] : computes recursively all CoulombH fac- 

tors for DSZ matrix Mat, perturbed to PMat, and any dimension vector strictly less 
than Nvcc; relies on the next two routines: 

• ListCoulombH[iVyec_,Q_] : returns returns a pair {ListH, ListC} where ListH is a list of 

CoulombH factors possibly appearing in the Poincare-Laurent polynomial Q of a quiver 
with dimension vector Nvcc, and ListC is the list of coefficients which multiply the 
monomials in Qsio^i^u) canonically associated to the H factors in Q. 

• SolveCoulombH[LjstH_,ListC_, R_ ] : returns a list of replacement rules for the CoulombH 

factors listed in ListH, by applying the minimal modification hypothesis to the coeffi- 
cients listed in ListC. The last argument is a replacement rule for CoulombH factors 
associated to subquivers. 

• MinimalModif[f_] : returns the symmetric Laurent polynomial which coincides with the 

Laurent expansion expansion of the symmetric rational function / at ?/ = 0, up to 
strictly positive powers of y. Here symmetric means invariant under y ^ 1/y. In 
practice, MinimalModif [i] evaluates the contour integral in [14], Eq 2.9 



by deforming the contour around into a sum of counters over all poles of f{u) and 
zeros of (1 — uy){l — u/y). This trick allows to compute (A.l) even if the order of the 
pole of f{y) at y = is unknown, which happens if SQuiverMultiplier is a formal 
variable. 

• SymplifyOmSbasis[f_]: replaces fls{l,y) 1 when 7 is a basis vector; 

• Symplif yOmSbasismult[f_] : replaces Qs{l,y) when 7 is a non-trivial multiple of a 

basis vector; 




(A.l) 
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• CoulombHNoLoopToZero[Mat_,f_]:sets to zero any H factor in f corresponding to sub- 

quivers without loop, assuming DSZ products aij = Mat[[i,j]] ; active only on 2-node 
subquivers if $QuiverTestLoop is set to False 

• OmTNoLoopToZero[Mat_,f_]:sets to zero any fitot factor in / corresponding to subquivers 

without loop, assuming DSZ products aij = Ka.t[[i,j]] ; active only on 2-node subquivers 
if SQuiverTestLoop is set to False 

• OmSNoLoopToZero[Mat_,f_]:sets to zero any Qs factor in f corresponding to subquivers 

without loop, assuming DSZ products a^j = Mat[[z, j]] ; active only on 2-node subquivers 
if SQuiverTestLoop is set to False 

• EvalCoulombH3[Mat_,f_ ] : evaluates any 3-center H factor with multiplicity vector {1, 1, 1} 

appearing in /. Not used in any routine so far. 

• DropFugacity[f_] : drops second argument in any Qsil,y) appearing in f, to be used at 

end of computation only; 

A. 5 Higgs branch formula 

• HiggsBranchFormula[Mat_,Cyec_,Arvec_,y_] : computes the Poincare-Laurent polynomial 

of a quiver with DSZ products a^j = Mat[[z, j]] (possibly rescaled by $QuiverMultiplier), 
dimension vector Nt = Nvec[[i]], FI parameters d = Cvec[[i]], using the Higgs branch 
formula (4.1). It is assumed, but not checked, that the quiver has no oriented loop; 

• StackInvariant[Afat_, Cvec_,iVyec_,j_ ]: gives the stack invariant (4.2) of a quiver with 

DSZ matrix ajj = Mat[[i,j]], possibly rescaled by an overall factor of $QuiverMultiplier, 
FI parameters Q = Cvec [[«]], dimension vector Ni = Nvec [[«]], using Reineke's formula 
(4.2); the answer is written in terms of unevaluated g-deformed factorials QFact[n,y]] 

• AbelianStackInvariant[Mat_,Cvec_,j_]: gives the Abelian stack invariant (4.7) of a 

quiver with DSZ matrix a^j = Mat[[i,j]], possibly rescaled by an overall factor of 
SQuiverMultiplier, FI parameters Q = Cvec[[i]], using Reineke's formula (4.2); co- 
incides with Stacklnvariant with Nvcc= {1, . . . 1} except that tests of marginal or 
threshold stability are performed (unless $QuiverVerbose is set to False); 

• QDef ormedFactorial[ji_,7_] :gives the g-deformed factorial [n,y]\ 

• EvalQFact[f_]: evaluates any QFact[n,y] appearing in f 
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A.6 Utilities 

• ListAllPartitionsCciarge yector_ ]: returns the list of unordered partitions {at} of the 

positive integer vector 7 as a sum of positive, non-zero integer vectors a^; 

• ListAllPartitionsMultCcharge yector_ ]: returns the list of unordered partitions {oj, mj} 

of the positive integer vector 7 as a sum of positive, non-zero integer vectors a, with 
multiplicity nii; 

• ListSubQuivers[Afvcc_]: gives a list of all dimension vectors less or equal to Nvec; 

• SubDSZ[Mat_, Cvec_,Li_ ]: gives the DSZ matrix of the subquiver made of vectors in list 

Li; 

• SyinmetryFactor[Li_] : gives the symmetry factor l/|Aut({ai, 0:2, ■ ■■ ,an}| for the list of 

charge vectors Li; 

• OmTRat[A''vcc_,y_ ] : gives the rational total invariant 0^1(7; Z/) in terms of Qtot{l',y)- 

Coincides with the latter if 7 is primitive. 

• OmTToOmS[f_ ]: expands out any Qtot{l]y) in / into H factors and fig's using(3.4); 

• OmToOmb[f_]: expresses any Q{'j;y) in / in terms of r2(7;?/)'s; 

• OmbToOin[f_ ]: expresses any ^(7; y) in / in terms of ^(7; ?/)'s; 

• HiggsGToOmb[iVvec_,y_]: Returns the (unevaluated) HN invariant GHiggs(7,l/) in terms of 

the rational refined indices ^(7;?/) using (4.5); 

• OmbToHiggsG[JVyec_, j_ ] : Returns the (unevaluated) rational refined index ^(7; y) in terms 

of the (unevaluated) stack invariants Giiiggs{l,y) using (4.1); 

• RandomCvec[Arvec_]: generates a random set of FI parameters d between -1 and 1, such 

that X;ONvec[[z]] = 0; 

• UnitStepWarn[A_] : gives 1 for a; > 0, for x < 0, and 1/2 if a; = 0. Produces a warning in 

this latter case, irrespective of the value of $QuiverVerbose. If so, the user is advised 
run the computation again with a different random perturbation. 
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